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ABSTRACT 



We use the recently calculated two-loop anomalous dimensions of current-current 
operators, QCD and electroweak penguin operators to construct the effective Hamil- 
tonian for AS = 1 transitions beyond the leading logarithmic approximation. We 
solve the renormalization group equations involving a s and a up to two-loop level 
and we give the numerical values of Wilson coefficient functions C(/x) beyond the 
leading logarithmic approximation in various renormalization schemes. Numerical re- 
sults for the Wilson coefficients in AB = 1 and AC = f Hamiltonians are also given. 
We discuss several aspects of renormalization scheme dependence and demonstrate 
the scheme independence of physical quantities. We stress that the scheme depen- 
dence of the Wilson coefficients C(/x) can only be cancelled by the one present in the 
hadronic matrix elements (Qi(fj,)). This requires also the calculation of 0(a) cor- 
rections to We propose a new semi-phenomenological approach to hadronic 
matrix elements which incorporates the data for CP-conserving K — > irir amplitudes 
and allows to determine the matrix elements of all (V — A) <g) (V — A) operators in 
any renormalization scheme. Our renormalization group analysis of all hadronic ma- 
trix elements (Qi{fJ>)} reveals certain interesting features. We compare critically our 
treatment of these matrix elements with those given in the literature. When matrix 
elements of dominant QCD penguin (Qe) and electroweak penguin (Qs) operators 
are kept fixed the effect of next-to-leading order corrections is to lower consider- 
ably e'/e in the 't Hooft-Veltman (HV) renormalization scheme with a smaller effect 
in the dimensional regularization scheme with anticommuting 75 (NDR). Taking 
m t = 130 GeV, Aj^g = 300 MeV and calculating (Qq) and (Qs) in the 1/N approach 
with m B (l GeV) = 175 MeV, we find in the NDR scheme e'/e = (6.7 ±2.6) x 10" 4 in 
agreement with the experimental findings of E731. We point out however that the 
increase of (Qq) by only a factor of two gives e'/e = (20.0 ± 6.5) x 10~ 4 in agreement 
with the result of NA31. The dependence of e'/e on Aj^t, m t and (Qq : s) is presented 
.A detailed anatomy of various contributions and comparison with the analyses of 
Rome and Dortmund groups are given. 
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1 Introduction 



The CP- violating ratio e'/e is governed by penguin contributions JT|, 0. Until 1989 there 
had been a general belief || f§ that e'/e is mainly described by the QCD-penguin contri- 
butions with the electroweak penguins (7- and Z°-penguins) playing only a secondary role. 
A new development in this field came in 1989 through the increase of m t and the work of 
Flynn and Randall || who pointed out that for large mt the electroweak penguin contri- 
butions could become important and, having opposite sign to QCD-penguin contributions, 
could considerably suppress e'/e. A detailed anatomy of e'/e in the presence of a heavy 
top quark has been done by Buchalla, Harlander and one of the authors 0. It has been 
found that with increasing m t the electroweak penguins could indeed compete with QCD 
penguins to cancel their contribution completely for m t ~ (9(220 GeV) so that e'/e ~ 
like in superweak theories. This rather surprising result has been confirmed in subsequent 
years by Paschos and Wu [0] and by Lusignoli, Maiani, Martinelli and Reina ||. 

On the experimental side, after heroic efforts on both sides of the Atlantic the situation 
is not yet conclusive || [Tcfl , 

„ , , . , f (23 ± 7) • KT 4 NA31 , . 

Re e'/e = I ) ' . . 1.1 

I (7.4 ±6.0) • 10- 4 E731 

While the result of NA31 clearly indicates a non-zero e'/e, the value of E731 is compat- 



ible with superweak theories in which e'/e = [|Tl|]. From the point of view of the analyses 
quoted above the NA31 result points toward a top quark mass of m t « (9(100 GeV) 
whereas the E731 result points toward m t ^ (9(150 GeV). Hopefully, in the next five years 
the experimental situation on e'/e will be clarified through the new measurements by the 
two collaborations at the 10~ 4 level and by experiments at the ^-factory in Frascati. 

Since the outcome of the fight between QCD and electroweak penguins is rather sensitive 
to various approximations used in refs. || [7], ||, it is important to improve the theoretical 
calculations of short distance (Wilson coefficient functions Cj(//)) and also the estimates 
of long distance (hadronic matrix elements (Qj (//)}) contributions. In the present paper 
we will make progress on both parts of the problem by calculating the Wilson coefficients 
beyond the leading logarithmic approximation and by extracting several hadronic matrix 
elements from the existing very accurate data on CP-conserving K — > tttt amplitudes. 
Needless to say that many results presented here are not only important for e'/e but can 
also be applied to all AS = 1, AC = 1 and AB = 1 decay processes. 

The present paper culminates our extensive studies of AF = 1 effective non-leptonic 
Hamiltonians beyond the leading logarithmic approximation |12|-||15||. In ref. fl2| two-loop 
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anomalous dimensions O(a^) of current- current operators have been calculated, confirming 



the previous results of Altarelli et al. jig] . This calculation has been extended to include 
QCD penguin operators and electroweak penguin operators in ref. |T5| , |TJj] where a two- 
loop 10 x 10 anomalous dimension matrix 0(a^) has been presented. Subsequently, the 
corresponding two-loop matrix of order 0(aa s ) has been calculated [|l^]. The latter is 
necessary for a consistent treatment of electroweak penguin operators beyond the leading 
logarithmic approximation. 

One of the aims of the present paper is to construct the effective Hamiltonian for 
AS = 1 transitions which incorporates the results of refs. [12|- [15]. This Hamiltonian will 
take the form 

/~i 10 C _ 

= QOi) = -|Q(//) T (7(//) , (1.2) 

where the Wilson coefficients Cj(/i) include leading and next-to-leading QCD corrections 
and leading order corrections in the electromagnetic coupling constant a. Here Qi <2 stand 
for current-current operators, Q3-6 for QCD penguin operators and Q7-10 for electroweak 
penguin operators. Explicit expressions for these operators will be given in section 2. A 
truncated effective Hamiltonian for Q1-6 with a = has been presented already by us in 



ref. [jT3 |, where an application to the AI = 1/2 rule and e'/e has been made. The present 



paper can be considered on one hand as a generalization of ref. |13|] to include electroweak 



penguin operators. On the other hand it can be considered as a generalization of the 
analyses |5], |], [7|, |J to include next-to-leading logarithmic effects. Whereas in refs. |5|, |6|, [7|, |[ 
the logarithms (a s t) n and at(a s t) n with t = ln(M^//i 2 ) have been summed, the present 
analysis includes also the logarithms a s (a s t) n and a(a s t) n . Moreover, we improve the 
treatment of hadronic matrix elements. 

The central questions which one would like to address in a paper like this one are as 
follows: 

• Do the complete next-to-leading corrections enhance or suppress the Wilson coeffi- 
cients of current-current operators ? 

• Do the complete next-to-leading corrections enhance or suppress the Wilson coeffi- 
cient functions of QCD penguin operators ? 

• Do these corrections enhance or suppress the Wilson coefficients of electroweak pen- 
guin operators relative to the ones of QCD penguin operators ? 

• What is the impact of all these effects on e'/e and the AI = 1/2 rule ? 
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How does e'/e depend on 



? 



Unfortunately, the answers to these questions are complicated by the fact that the 
coefficients Cj(/i) depend at the next-to- leading level on the renormalization scheme for the 
operators Qi, and only after the hadronic matrix elements (Qi(/j,)) have been calculated at 
the same level, can this scheme dependence be removed. This requires also the calculation 
of 0(a) corrections to (QiQjL)). Thus the size of next-to-leading order corrections to Cj(/x) 
calculated in the NDR scheme (anticommuting 75 in D 7^ 4 space-time dimensions) differs 
from the corresponding corrections found in the HV scheme (non-anticommuting 75 in D ^ 
4 space-time dimensions) . In order to cancel these dependences one would have to calculate 
(Qi(fj)) in a method sensitive to such a scheme dependence. To our knowledge only QCD 
sum rules could be useful in this respect at present but future lattice calculations also could 
overcome this difficulty, at least in principle. In spite of this our short distance calculation 
provides an important and necessary step towards a more accurate estimate of weak decay 
amplitudes than the one given merely by the leading logarithmic approximation. Moreover, 
in order to consistently study the m t -dependence of weak amplitudes it is mandatory to 
go beyond the leading logarithmic approximation as will be evident from what follows. 

The scheme dependence just discussed is an important point which we would like to 
address here. To this end we will evaluate the Wilson coefficient functions in three different 
schemes and we will investigate how the scheme dependence affects the answers to the set 
of questions stated above if no proper scheme dependence is incorporated in the matrix 
elements (Qi(fi)). In particular, we find as a rather unexpected result that the next- 
to-leading QCD corrections to Wilson coefficients of current-current operators Q\ and Q2 
calculated in the renormalization schemes considered in fact suppress slightly the AI = 1/2 



transitions and enhance the AI = ~ transitions. We explain why the authors of refs. |12 
and |16| who performed correct calculations reached conclusions opposite to ours. 



Other not unrelated issues are the /^-dependence of (Qi(p)) and the actual values of 
these matrix elements. For this reason we have made a critical comparison of the exist- 
ing estimates of (Qi(fi)) based on the vacuum insertion, the 1/iV-expansion and lattice 
methods. In this connection, we point out that the evolution of (Qi(fi)) with \i can be 
calculated within renormalization group improved perturbation theory. Performing to our 
knowledge the first detailed renormalization group analysis of (Qi(f-i)) we find that due 
to the asymmetry of the anomalous dimension matrices the /i-dependences of (Qi(f^)) are 
quite different from the ones of Cj(/x). This allowed us to answer the following question: 

• Should the parameters Bi entering various (Qi(n)) be /i-dependent in order to be 
consistent with the true QCD evolution ? 
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Rather unexpectedly, we have found that the Bi parameters of all (V — A) <g> (V + A) 
operators Q 5 - Qg show only very weak /x-dependence for /j, ^ 1 GeV supporting to some 
extent the structure of these matrix elements which is common to the three approaches in 
question. On the other hand the Bi parameters of (V — A)®(V — A) operators show sizeable 
/x-dependence which is fully missing in the vacuum insertion approach, partly present in 
the 1/N approach and not yet calculable in the lattice approach. 

In view of the latter problem and the fact that the present methods are not yet capable 
of addressing properly the renormalization scheme dependence of (Qi(fj)), nor fully explain 
the known data on CP-conserving K — > tttt decays, we have developed a phenomenological 
framework in which the matrix elements (Qi(/j)) necessary for the calculation of e'/e are 
extracted as far as possible from the data on CP-conserving K — > tttt decays. In this 
framework the central role is played by the calculable scheme and /z-dependent coefficients 
Ci(n) and the scheme and fx- independence of measured K — > tttt amplitudes gives us 
automatically the scheme and /i-dependence of (Qi(fi)) which are in addition consistent 
with the CP-conserving K — > tttt data. Thus the AI =1/2 rule is incorporated in our 
approach from the beginning. We verify that the scheme and /x-dependences of {Qi{n)) 
found this way are to a very good approximation consistent with the ones given by the 
renormalization group. 

Our approach has the four basic parameters 

Ams, P4 V2) (m c ), 4 1/2) (m c ), Pf 2) (m c ). (1.3) 

The choice /J, — m c turns out to be very convenient but is not necessary. The last three pa- 
rameterize hadronic matrix elements for K — > tttt, and will be defined in section 5.2. 
B^ 2 \m c ) parameterizes all matrix elements (<5i)o of (V — A) <g> (V — A) operators. 
B^ 2 \m c ) and B^ 2 \m c ) parameterize the (V — A) <g) (V + A) operators. In our phe- 
nomenological approach all matrix elements are given in terms of these four parameters in 
such a way that for each set of these parameters all data on CP-conserving K — > tttt am- 
plitudes are reproduced. Moreover, all the matrix elements (Qj(A i ))2 of (V — A) <g> (V — A) 
operators (i = 1, 2, 9, 10) are given in our approach as numerical functions of only A^, \x 
and the renormalization scheme considered. With all these matrix elements at hand we 
can predict the CP- violating quantity 

ImA t 



e' 



= 10 

e 

where A t = V* s V tA and 



1.7 • 10- 4 



p(l/2) _ p(3/2) 



;i.4) 



= 4 1/2) + 4 1/2) Bi 1/2) + 4 1/2) Bi 1/2) , (1.5) 
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pirn = 4 3/2) + 4 3/2) Bf /2) 



(1.6) 



We calculate the coefficients aj as functions of A-^ and m t for the leading order and the 
two renormalization schemes (NDR, HV) considered. 

We hope very much that this novel anatomy of various aspects of e'/e including next- 
to-leading order corrections and an improved treatment of hadronic matrix elements will 
shed some light on the present uncertainties and on the theoretical prediction of e' / e in the 
Standard Model (SM). By doing this we will attempt to clarify the discrepancy between 
recent theoretical estimates of e'/e which bears some similarity in it to the discrepancy 
concerning experimental results on e'/e given in eq. (1J-). Whereas the Dortmund group 
|T7| , |T8"| found last year e'/e to be close to the NA31 result, the very recent analysis of the 
Rome group [|1£]] favours the experimental result of E731. 

Our paper is organized as follows: In section 2, we will present general formulae for 
Ci(fi) beyond the leading logarithmic approximation. These formulae generalize the ones 
given in refs. [|li| and ||. In section 3, the renormalization scheme dependence will be 
discussed. In section 4, we calculate the coefficients Cj(Mw) and use the known one-loop 
0, [2~0]-[2~7] and two-loop anomalous dimension matrices [|T2]]-|]n| to find Cj(/i) in various 



renormalization schemes. In section 5, we collect and compare various existing theoret- 
ical estimates of hadronic matrix elements. In section 6, we present the details of our 
phenomenological approach to the matrix elements, we calculate the matrix elements as 
functions of the parameters in eq. ( |1.3|) and compare the results with the ones obtained in 
section 5. Here, we also present the renormalization group evolution of (Qi(fi)) advertised 
earlier in this paper. In section 7, we recall the CKM-parameters necessary for our nu- 
merical analysis of e'/e and we give the AS = 2 Hamiltonian which we used to determine 
the CKM-phase 5 from the experimental value of the parameter e describing indirect CP- 
violation in the .fT-system. In section 8, we give the basic expressions for e'/e in terms of 
Ci(fi) and {Qi(fi)), and we address all the issues mentioned above. In particular we give 
the results for e'/e as a function of Aj^g, m t and the matrix elements (Qe)o and (Qs}2- 
For completeness, we give in section 9, the numerical values for Ci(fi) relevant for AB = 1 
transitions. We end our paper with a list of main messages coming from our new anatomy 
of e'/e. 
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2 Basic Formulae for Wilson Coefficient Functions 



2.1 Operators 

In order to illustrate the scheme dependence of Wilson coefficient functions we will use two 
bases of operators. Our main basis is given as follows 



Qx = 


(s a up) y _ A (u p d a ) v _ A , 


Q2 = 


(su) v _ A (ud) v _ A , 


Q 3 = 


(* d )v-A (QQ)v-A ' 
q 


Q 4 = 


{s a d l 3) v _ A ^2(q l3 q a ) v _ A , 
q 


Q 5 = 


(sd) v _ A (qq)v+A » 


Q 6 = 


{s a d(3) v _ A ^2{q(3q a ) v+A , 
q 


Qi = 


3 

(^)v-A £ C 1 (W)v+A > 
Z <? 


Qs = 


3 

(s a d ) y _ A ^ e g (q q a ) y+A 
z q 


Q9 = 


3 

(sd\^ A Y,e q {qq\_ A , 


Qio = 


3 

n (M/Ov-A X! e 9 (Wcdy-A 
z g 



(2.1) 



where a, /3 denote colour indices (a, /3 = l,...,iV) and e g are quark charges. We omit 
the colour indices for the colour singlet operators. (V ± A) refer to 7 M (1 ± 75). This basis 
closes under QCD and QED renormalization. 

For n < m c the sums over quark flavours in (|2.1| ) run over u, d and s. However, when 
m b > fi > m c is considered also q = c has to be included. Moreover, in this case two 
additional current-current operators have to be taken into account, 

QI = (s a c /3 ) y _ A {cpd a ) y _ A , QI = (sc) v _ A (ccf) v _ A . (2.2) 

It should be stressed that this basis differs from the one used by Gilman and Wise [|23| 
and by the Rome group [|19[] in the form for Qx, Q 2 , QI and Q\ where the corresponding 
Fierz conjugates have been adopted 

Qx = (sd) y _ A (■uw)v-a > = (s a d(3) Y _ A (uf3U a ) Y _ A , (2.3) 
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Ql = (sd) v _ A (cc) v _ A , Q c 2 = {s a d p ) y _ A (c /3 c a ) v _ A . (2.4) 

We will refer to the basis of refs. [19| as the second basis. Both are equally good although 
we somewhat prefer the basis (|2.1|) , because there Q2 is taken in the colour singlet form as 
it appears in the tree level Hamiltonian. 

Now, as we have stressed in refs. [ 13| , |i~4 , pTSfl the two-loop anomalous dimension matrices 



in the two bases differ when calculated in the NDR scheme. Consequently, the correspond- 
ing coefficients Cj(/i) differ. We will refer to these two possibilities as schemes NDR and 
NDR, respectively. On the other hand in the HV scheme it is immaterial whether the first 
or the second basis is used but the corresponding HV Wilson coefficient functions differ 



from those obtained in the NDR and NDR schemes. 

In what follows we will present general formulae for Wilson coefficients valid in any 
scheme leaving the discussion of scheme dependence to section 3. 

We should also remark that we do not include in our analysis the operators 

Q11 = (^)v-a (») V _ A , Q12 = (S6) V _ A (M) v _ A , (2.5) 

which have been taken into account in ref. || [7]]. Even at leading order their effect on e'/e 
is smaller than 1% and in order not to complicate the analysis further, we have decided to 
drop them. This was also the strategy of the Rome group ]TP| . 



2.2 Renormalization Group Equations 

The renormalization group equation for C(/z) is given by 



9 n , X 9 



C(^f,S 2 , a) = 7 V, <*) <?(#, 9 2 , <*) , (2.6) 



where (3{g) is the QCD beta function 

M = ~ ^ ~ Pi-^ - Pu^p , (2-7) 

with 

/3o = H-|/, A = 102-|/, = -!(„+!), (2.8) 

and / = u + d denoting the number of active flavours, u and d being the number of w-type 
and d-type flavours respectively. To the order considered in this paper, in principle we 
also have to take into account the third term in eq. [2.8| . However, we have checked that 
its contribution to the coefficient functions is negligible, and we will drop it for the rest of 
this work. 
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In what follows, we also neglect the running of the electromagnetic coupling constant 
a. This is a very good approximation because only scales 1 GeV < \i < are involved 
in our analysis. In the numerical analysis we will take a = 1/128. For the effective QCD 
coupling constant we will use 



An 



1 - 



(51 HQi/A}) 



(2.9) 



with Af in the MS scheme. Demanding the continuity of af (Q) at quark thresholds in 
the form 

af (m c ) = a W (m c ) , af (m b ) = af (m b ) , (2.10) 

gives relations between the various values for Af. In what follows we will denote Ayg? = 
A4 = Aqcd- 

Next, j(g 2 , a) is the full 10 x 10 anomalous dimension matrix which we expand in the 
following way 

' " ' f2.ll) 



7(^ 2 ,a) = 7s(/) + ^r(^) 



where 



and 
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47T 



(47T) 



2 7 S 5 



= 7 f + ^ 



(2.12) 



(2.13) 



Explicit expressions for (7f ,7f ) and {yjP\jj$) can be found in refs. |L4| and |f[5| , 
respectively. They will not be repeated here. 

The solution of the renormalization group equation ( |2.6| ) is given by 



C{^,g ,a 



T 9 exp / d# 



(2.14) 



where T 9 denotes the ordering in the QCD coupling constant such that the couplings 
increase from right to left. 



2.3 Generalized Evolution Matrix 

Let us next introduce a generalized evolution matrix from m<i down to m\ < 

rg(mi) 7 T (g' 2 , a) 
U(m 1 ,m 2 ,a) = T s exp / d# — — — — , 

Jg(m 2 ) p[g ) 



(2.15) 



S 



which is the generalization of eq. (2.3) of ref. ||13|| . Eq. ( |2.14j ) can be now in a short hand 
notation written as 

d(jj) = U(fi, M w , a)(7(Mw) • (2.16) 
The matrix U(m 1 ,m 2 ,a) can be decomposed as follows 

(2.17) 



a: 



U(m ll m 2 ,a) = U(m 1 ,m 2 ) + —R(m 1 ,m 2 ) , 

Alt 

where t/(m 1 ,m 2 ) describes pure QCD evolution and R(mi,m 2 ) the additional evolution in 
the presence of the electromagnetic interaction. U(mi, m 2 ) sums the logarithms (a s t) n and 
a s (a s t) n with t = \n(m\/m\) whereas i?(m l5 m 2 ) sums the logarithms t(a s t) n and (a s t) n . 
f/(mi,m 2 ) has been given in ref. The leading order formula for R(mi,m 2 ) can be 

found in ref. || except that here we used a different overall normalization (an additional 
factor —An in R). For completeness, we will now recall the formula for f/(m 1 ,m 2 ) and 
subsequently, we will derive a corresponding expression for R(mi,m 2 ). 



2.4 The Pure QCD Evolution Matrix U(m 1 ,m 2 ) 

As shown in [[n| the pure QCD evolution can be written as 

A a 

U(m 1 ,m 2 ) — H — 



(2.18) 



An j V 47r , 

where U^(m%, m 2 ) denotes the evolution matrix in the leading logarithmic approximation 
and J summarizes the next-to-leading correction to this evolution. If 

(7 s (0) )d = V-^V , G = V-'t )T V , (2.19) 

where (j^)d denotes a diagonal matrix whose diagonal elements are the components of 
the vector then 



f> (0) (mi,m 2 ) = V 
For the matrix J one gets 



aJm 2 



'Y 



v- 1 



D 



a s (mi) / 



with 



(0) 



2/3 



where the elements of S are given by 



Si 



Sij 7. 



(o) Pi 



G; 



2(3 + ^-^ 



(2.20) 



(2.21) 



(2.22) 



with 7^ denoting the components of and the elements of G. Although eq. ( |2.22|) 
can develop singularities for certain combinations of the %}, the physically relevant evo- 
lution matrix ( |2.18j ) always remains finite after proper combination of relevant terms. 
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2.5 The Evolution Matrix R(mi, 777,2) 

Inserting ( [2.1 1|) into ( 2.15|) , we find a general expression for R(mi,m 2 ) 



R(m 1} m 2 ) 



g(mi) U(m 1 ,m')T T (g')U(m',m 2 ) 

d-9 7TTT\ 

9{m 2 ) /3{g') 



(2.23) 



where f/(m 1 ,m 2 ) is given by ( |2.18| ), T(g) has been defined in ( |2.13| ) and g' = g(m'). 

It is instructive to discuss the leading order in fl2.23|) which is obtained by keeping only 
the first terms in (|2.7|) and (|2.13|) and setting J = in U (mi, 7712) of eq. (|2.18|) . We then 
find 

R {0) ( mi ,m 2 ) = V k {0 \m 1 ,m 2 )V' 1 , (2.24) 
Po 



where the matrix i^^(m 1 ,m 2 ) is given by 



(k^\ mil m 2 )) l3 = M\f 



a s (mi) 



dal 



a s (m. 2 ) 



^ \ a ' J_ ( asijn 2 ) 
a s (mi)/ 



(2.25) 



with aj being the components of the vector a defined in eq. (p. 20 ). Moreover, the matrix 
is given by 

= V~ 1 ^ T V. (2.26) 



Straightforward integration gives 



M. 



(o) 



'j 



(m 2 ) \ a3 



a s (mi)J a s (mi) U s (mi)l a s (m 2 ) 



^2.27) 



which is precisely eq. (3.23) of ref. |J. 

Similarly to eq. (|2.22|) , we note an apparent singularity in the element (7, 8) for which 
a-j = a$ + 1 when / = 3. However, the expression in parenthesis also vanishes in this 
case and no singularity is present. For numerical calculations it is safer to use in the case 



a,- = a. 



+ 1 the following formula obtained directly from eq.( [2.25|) 



(0) 



13 a B (mi) Va s (mi) 



a s (m 2 )\ aj aig (mi) 



a s (m 2 ) 



(2.28) 



It is straightforward to generalize ( |2.24j ) beyond the leading logarithmic approximation. 
We find 



.R(mi, m 2 ) 



2ty 



— V k{m x ,m 2 )V~ x = #(mi,m 2 ) + # (1) (mi, m 2 ) , 
Po 



where 



1 



k(mi,m 2 ) = k^\m 1 ,m 2 ) + — V-K~j (mi,m 2 ) . 

47r i=i 



(2.29) 



^2.30) 
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The next-to- leading corrections originating in J 7^ 0, 7^ and (3\ 7^ are repre 

sented by the matrices i^f 1 ' ) (m 1 ,m 2 ) and more globally by (m 1? m 2 ). 
Let us introduce 

f(i) 



.(1)T _ A *(0)T 
/se /o /e ' 



and 



M(i) = y-i (f«+[7f T ,j])v. 



The matrices X,- (7711,7712) are then given as follows 



(i) 



m 1 ,m 2 , 



\a s (mi)/ Us( m i)/ J ' J 
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\a B (mi)J a s (m2) 



(2.31) 
(2.32) 



(2.33) 



k!p(mi, m 2 ) 



a s (m 2 ) k(°\mi,m 2 ) S , 
a s (mi) S (mi, ra 2 ) • 



(2.34) 
(2.35) 



Comparison of ( |2.27| ) with ( |2.33| ) and the additional factors of a s (rrii) in ( |2.34 ) and 
Q2.35| ) make it clear that BP-'(mi, m 2 ) is by one logarithm lower or by one order in a s 
higher than J?^(m 1 ,m 2 ). For later purposes it will be useful to introduce "a s -counting" 
in which i^ ) (mi, 777,2) is 0(l/a s ) and R^(mi,m 2 ) is 0(1). 



2.6 A Different Form of the Evolution Matrix 

The evolution matrix U (mi, 777.2, <x) can also be written in a form which nicely generalizes 
the pure QCD evolution matrix of eq. ( |2.18| ), 



where 



W(rm) 
W'(m 2 ) 



U(m 1 ,m 2 ,a) = W(mi)U^\m 1 ,m 2 )W\m 2 ) 

a B (mi) 



-J 

a 



1 + 



-J e )(l 



J, 



a 



An 1 ' a s (mi 
a s (m 2 ) a 



47T 



(2.36) 

(2.37) 
(2.38) 



0:47712 J 

This is in fact the form used in ref. Ul~9| , where only implicit equations for the matrices 
J s , J e , and J se have been given. Here, we give explicit expressions for the Jj which can be 



easily found from the previous subsections. J s is simply given by J of ( |2.21|) and J e , J se 
are found to be 

J e = V S e V' 1 , J se = V S se V~ l , (2.39) 
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where 



I M, 



(o) 

ij 



{Se)t > 2(3 (1 + aj - a, 



1 



(S se )ij — n r> 



M, 



(i) 



2(3 (a 7 - - a { 



(2.40) 



with M$ and M^ 1} defined in eqs. ( |2.26| ) and ( |2.32| ), respectively. We note that S e can 



develop singularities for certain combinations of aj and a«, and S se is singular for i = j. 
All these singularities cancel in the final expression for If (mi, rri2, a) when all contributing 
terms are combined. This is seen for instance in (|2.33| ). Thus although formula ( |2.36| ) 
is more elegant than the evolution matrix presented in sections 2.4 and 2.5, the former 
formulation in terms of the evolution matrices U(mi,m2) and R(mi,rri2), which is non- 
singular at all stages, is more suitable for numerical calculations. 



2.7 The Initial Conditions C(M W ) 

In order to complete the analysis, we have to discuss the structure of C(Mw). This 
discussion generalizes the one given in section 2 of ref. fll3| . In order to find C(Mw) the 
one-loop current- current and penguin diagrams of fig. [1] with the full W and Z propagators 
and internal top quark exchanges have to be calculated first. Subsequently, the result of 
this calculation should be expressed in terms of matrix elements (Q(M^)). The latter are 

— # 

found by inserting the operators Q in the one-loop current-current and penguin diagrams of 
fig. H and calculating the finite contributions in some renormalization scheme. We note the 
absence of Z°-contributions to (Q(M^)). Indeed, in the mass independent renormalization 
schemes used here such contributions are zero. 
If 



T=(QW T ) 



f (o) + 



a s (Mi 



s 4tt e 



47T 



(g(M w ) T )C(M 1 



w 



(2.41) 



denotes the result of evaluating the diagrams of fig. [l| with (Q^) being the tree level matrix 
elements and 

V o) > 



obtained from fig. 0, 
C(M W ) 



(Q(M W )) 
then 

a s (M w 



~ a s (M w 
1 H : r 



47T 



s + -r- r e 



47T 



f (0) + 



An 



a 

+ — 

47T 



(2.42) 



^2.43) 



It should be stressed that r s and f e depend on the assumptions made about 

the properties of the external lines in figs. 1 and 2, i.e., on the infrared structure of the 
theory. This dependence cancels however in (|2.43|) so that 0(a s ) and 0(a) corrections 
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in ( |2.43j) do not depend on external states as it should be. They depend however on the 
renormalization scheme through the matrices f s and f e . 

The formulae ( |2.41| )-( p.43| ) are general and apply also to situations in which several 
operators are present in the limit a = a s = represented by the first term in ( |2.43| ). In 
the case at hand only the operator Q 2 is present in this limit and the vector is the 
transposed of (0, 1, 0, . . . , 0). This implies that only ((^(-^w)) has to be calculated in order 
to find C(Mw). This is evident from (|2.43|) and the discussion in appendix C of ref. [p~3 



Yet, this more general formulation is very useful for the analysis of scheme dependences as 
we will see below. 



2.8 Final Formulae for Wilson Coefficient Functions 

Let us define 



a 



cm = cm + ^ cm 



(2.44) 



Inserting (|2~T7| ) and ( glj) into (gig ), we find 



cm = (i + u^, m w ) ( ?v + « - <M^l { rJ + j)f^) , 

\ 4lT J \ 47T 47T I 



(2.45) 



and 



^•o*) = f> {o) (^Mw)[f«-ef(°) 

+ R (1 \fi,M w )f 



f {0 ) + «s(M w ) 



47T 



?(0) 



(2.46) 



3 Renormalization Scheme Dependence 



3.1 General Consistency Relations 

The two-loop anomalous dimensions depend on the renormalization scheme for operators 



we 



and in particular on the treatment of 75 in D 7^ 4 dimensions. In refs. [|T3| and |T5 
have derived two relations between 7W and 7^ calculated in two different renormalization 
schemes a and b. 



a (1) ) 6 



(7 s (1) ) a + 



Af s , 7s (0) +2A,Af s 



{1 \ 

se Jo 



(7. 



se J a 



+ 



Af s ,7 



(0) 



Ar e ,7 



(0) 



(3.1) 
(3.2) 
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with 



An 



and fi denned in ( p. 42 ). It follows that the combinations 



7 



(i) 



(0) 



2(3 r s 



(3.3) 



(3.4) 



and 



/ se 



' s, Je 



(0) 



(3.5) 



are independent of the renormalization scheme considered. We should also point out that 
Af s and Af e do not depend on the infrared structure of the theory, whereas this is not the 
case for f s and has been stressed in section 2. 

Relations (p. l|)-(pT4|) turned out to be very useful in testing the two-loop anomalous 
dimension matrices found in refs. and [15|. They also play a central role in demon- 
strating the scheme independence of physical quantities. For the pure QCD case this has 
been discussed in ref. ||13||. Below we will extend these considerations to 0(a) corrections. 



3.2 Relations Between Coefficient Functions 



The relations between coefficient functions obtained in two different schemes a and b are 
entirely given in terms of the matrices Afj of eq. (|3.3|) which follow from the evaluation of 
the one-loop diagrams of fig. |2|. 
One has 



CM 



47T 



Af 



a 



T-Af e T 

47T 



with the corresponding relation for the hadronic matrix elements 

a s (jd) 



{Q{p))i = (Q0*))i 



i + 



- — Aff + ^AfJ 

47T 47T 



(3.6) 



(3.7) 



The scheme independence of the effective Hamiltonian ( |1.2| ) follows directly from these 
relations. 

The matrices Af s and Af e relating HV and NDR schemes can be found in sections 3.3 
and 3.4 of ref. [14| and in section 6.2 of ref. [0, respectively. The corresponding matrices 
relating NDR and NDR schemes are given in eq. (5.14) of ref. [JTSj . See however eq. ( f4.25|) 
of the present paper. 
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3.3 Scheme Independence of Physical Quantities 



Let us return to the general formulae (|2.45 ) and ( |2.46| ) for the coefficients C(fi). In ref. [13 



we have shown that the combination f J + J in eq. Q2.45| ) is renormalization scheme in- 
dependent, implying that the next-to-leading QCD corrections to C s (n) are only scheme 
dependent through the matrix J at the lower end of the evolution which involves a s (/j,). 
This scheme dependence is cancelled by the one of 0(a s ) corrections to (Q(fJL)) so that 
the resulting decay amplitudes do not depend on the renormalization scheme as it should 
be. We would like now to demonstrate that an analogous proof can be made in the case 
of C e (fj) given by ([2.46 ) . This can certainly be done, but it turns out that a much more 
elegant proof can be made by using the form of ( |2.36fl . 

First we want to demonstrate the scheme independence of terms at the upper end of 
the evolution. Inserting (|2.36| ) and ( 2.43 ) into (|2.16|) and recalling that T^°\ and J e 
are scheme independent we find that it is sufficient to consider the combination 



X 




a r 

47T 



T 



(3.8) 



The first part of this expression is scheme independent due to scheme independence of J c 
and rj + J s as shown in |13||. It remains therefore to show that also + J se is scheme 
independent. 

Now the scheme independence of W se and + J s allows to extract the scheme dependent 
part of M (1) : 

Mjp (scheme dependent) = (7$ - 7$) (P e )y ( 3 - 9 ) 



where P c = V 1 fj V. Therefore using (|2.39| ) and ( 2.40 ) we conclude that the scheme 



dependent part of J se is simply — rj which ends the proof. 

In a very similar way one can show that the scheme dependence around the \x scale 
present in W(fi) of (p . 3 7|) is cancelled by the scheme dependence of the 0{a s {jj)) and 
0(a) corrections to the matrix elements (Q T (fi)). This ends the proof of the scheme 
independence of physical amplitudes. 



4 AS = 1 Hamiltonian Beyond Leading Logarithms 
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4.1 General Remarks 



In the absence of strong and electromagnetic interactions the effective tree-level Hamilto- 
nian for AS = 1 non-leptonic decays can be written as follows 

H eS (AS = 1) = ^§KdK* s {(1 - r)(Q u 2 - Ql) + t(Q u 2 - Ql)} , (4.1) 

where 

r = - I and Q q 2 = {sq) V - A {qd)v-Ai q = u,c,t. (4.2) 



If the NDR scheme is adopted the Fierz conjugates Q\ of ( F2TB| ) and ( |2.4| ) are present in 
( |4.1|) . For the HV scheme it is immaterial whether Q\ or Q\ are used. 

The inclusion of QCD corrections and the evolution of 7i e s from Mw down to « Mw 
has been discussed in the leading logarithmic approximation at length in |3|. The next- 
to- leading order formulae involving Qi - Q e have already been given in ref. [HJ. Here, 



we generalize the results of |L3[ to include also the electroweak penguin operators Q7 - 
Qio together with 0(a) corrections using the formalism developed in previous sections. A 
similar though less explicit generalization has been made by the authors of ref. |19|] . In the 
latter paper however mainly the Wilson coefficients relevant for e' / e have been considered. 
Here, we present also the results for the coefficients relevant for the real parts of K — > ttti 
amplitudes which enter the discussion of the AI =1/2 rule. Moreover, it will be useful 
to discuss the cases /x < m c and m c < \i < separately, as well as the interesting case 
fi = m c on the border line of effective 4- and 3-flavour theories. This case is simultaneously 
relevant for AC = 1 transitions after proper replacement of quark flavours and is a useful 
renormalization point for our approach to hadronic matrix elements. The case \x ~ 0{m^), 
relevant for AB = 1 decays, is considered in section 9. 

4.2 Master Formulae for Wilson Coefficient Functions (// < m c ) 



The transformation of 7i e s of fl4.1|) from Mw down to /x < m c involves 5-, 4- and 3-flavour 
effective theories and can be compactly described as follows 

10 10 

Ql - Q2 - E ^0) Qi, Ql- Qi - E Q* » ( 4 -3) 

i=l i=l 



where Qi are the operators given in (2.1), with q = u,d,s. Consequently, we find the 
familiar expression 
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H cS (AS = 1) = %KdK*s E + r Viiv)) Qi(v) , (4-4) 

v2 i=1 
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where 

Vi{y) = Vi(fx) - Zi{n) • (4.5) 

The coefficients Zi and V{ are the components of the ten dimensional column vectors z 
and v 

z{n) = U 3 (fx, m c , a)z(m c ) , (4.6) 

and 

v(li) = U 3 (fi, m c , a)M(m c )U 4 (m c , m b , a)M(m h )U 5 (m h , M w , a)C(M w ) ■ (4.7) 

The full evolution matrices £7/ (mi, rri2, a) are given by ( |2.17| ) with / denoting the number 
of effective flavours. M(rrii) is the matching matrix at quark threshold mj. Its origin and 
structure is discussed separately in section 4.3. 

Using fl4.6p and ( |4.7| ), care must be taken to remove consistently higher order terms 
in a s and a resulting from the multiplication of the evolution matrices and initial values. 
To this end one should remember our discussion of the "a s -counting" in the case of the 
evolution matrix i?(m 1 ,m 2 ). From this counting it is evident that e.g. terms involving 
products JRS l > should be discarded to the order considered here. 

The initial values C(Mw) necessary for the evaluation of Vi(fi) are found according to 
the procedure outlined in section 2. For the case a = they have already been given in 
111 311. For the NDR scheme we find 



Ci(M w 
C 2 (M W 
C 3 (M W 

C 4 (M W 
C 5 (M W 

C 6 (M W 

CV(M W 
C 8 (M W 
C 9 (M W 
Cio(M w 



1 



47T 

| Qs(Mw) 

4:71 

a s (M w ) 



gNDR , a £?NDR 
47T '''' 



D NDR , " R NDR 



24tt 
a s (M w ) 
8tt 

a s (M w 



24tt 
a s (M w ) 



■E{x t ) 



a 

47r" 
a; 1 

67r sin 2 6 1 



[2B{x t )+C{x t )\ 



57T 



a 

Qtt 
0, 

a 

6tt 
0. 



4C{x t )+D{x t ) 
AC{x t ) + D(x t ) 



sin 2 W 



'10B{x t )-4C{x t )) 



(4.8) 
(4.9) 
(4.10) 
(4.11) 
(4.12) 

(4.13) 

(4.14) 
(4.15) 
(4.16) 
(4.17) 
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where 



rNDR 
e,l 

E(x t ) 



11/2 
0, 



E(x t )--, 



rNDR 

rNDR 

-°e,2 



-11/6, 
-35/18. 



D(x t ) = D(x t ) - - , 



m t 



and 



C(x) 
D(x) 
E(x) 



x 



X 

6 



+ 



In a; 



— Inx + 
9 



— In x + 

3 



a; In a; 

yX ~ l) 2 

3x + 2 
H (x - l) 2 

-19x 3 + 25x 2 , x 2 (5x 2 -2x-6) 

36(x- l) 3 " 
x(18 — llx — x 
12(1 -x) 3 



18(x- l) 4 



Inx , 



+ 



x 2 (15 - 16x + 4x 2 
6(1 -s) 4 



Inx . 



(4.18) 

(4.19) 
(4.20) 

(4.21) 
(4.22) 
(4.23) 
(4.24) 



Comparing with eqs. (3.26)-(3.35) of ref. 0, we observe that the initial conditions for 
the penguin operators given here differ from the ones used in ref. 0. There the leading 
order evolution matrices have been used but in the initial conditions the m t -dependent 
next-to-leading terms have been already taken into account. These m t -dependent terms 
represented by the functions B(x t ), C(x t ), D(x t ) and E(x t ) are renormalization scheme 
independent. Yet, the initial conditions may also contain scheme dependent next-to-leading 
terms which of course were not present in the analysis of ref. . In the NDR scheme these 
scheme dependent terms enter the initial conditions for Qi and Q2 and also modify the 
functions E(xt) and D(xt) by constant terms as given in ( |4.19| ). 

The initial values for the NDR and HV schemes can be found by using the relations 
( ^.6[ ). It is then interesting to observe that in these schemes the additional terms in ( f4.19p 
are absent and the initial conditions for Qi and Q2 are changed. The corresponding two- 
loop anomalous dimensions also change so that the final result for Ti, e g remains scheme 
independent. This clarifies the observation made by Flynn and Randall || that the initial 
conditions for penguin operators do depend on the form of the operators. As already em- 
phasized in |13| , the inclusion of two-loop anomalous dimensions cancels this dependence. 
Compared to ref. |L3|, we have however slightly changed our treatment of the HV scheme. 
Since in the HV scheme the anomalous dimension of the weak current 7j is non-zero at 
O{ol) [0, an additional contribution arises which can be either included in the anomalous 
dimension matrix as 7 S — 27j 1 \Y1\ or in the initial condition. In |TB|] we choose to work 
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with the latter treatment, which however introduces an unnecessarily large scheme depen- 
dence at the lower end of the evolution. Therefore, here we have included the anomalous 
dimension of the weak current in the anomalous dimension matrix. This also changes the 



initial conditions in the HV scheme, as was already discussed in ref. [13]. Moreover, Ar. 



given in [|13|, [14| has to be changed to 

Af s -> Af s - 4 C F 1 . (4.25) 

In order to calculate one has to consider the difference Q 2 — Q 2 . Then due to 

the GIM mechanism the coefficients Zi(fi) of penguin operators Qi, i ^ 1,2 are zero in 
5- and 4-flavour theories. The evolution in this range of // (// > m c ) involves then only 
current-current operators = Q\i and Q\ 2 with the initial conditions 

, M v _ «s(M w ) NDR a(M w ) rND r 9R n 
Zi(M w J - — — — B sl H — — B el , (4.2b ) 

~ (M \ - 11 ft s(M W ) RNDR , a ( M w) R NDR (a 97 N 

z 2 (M w J - H ^j- — i3 s , 2 H ^" — ^ e ,2 > ( 4 -27j 



as in ( |4.8| ) and (|4.9). 2 and Qi 2 do not mix under renormalization with each other and 
the initial conditions for the coefficients of Q\ 2 to be denoted by zl 2 (fi) are identical to 
the ones given above. We then find 

Zl{mc l ) = U,(m c ,m h ,a) M(m h ) U 5 (m b ,M w ,a) ( Zl{ ^\ ) , (4.28) 
Z2{m c ) J y z 2 {M w ) J 

where this time the evolution matrices contain only the (2, 2) anomalous dimension 
matrices describing the mixing between current- current operators. When the charm quark 
is integrated out the operators QI 2 disappear from the effective Hamiltonian and the 
coefficients Zi(fi), i ^ 1,2 for penguin operators become non-zero. In order to calculate 
Zi(m c ) for penguin operators a proper matching between effective 4- and 3-quark theories 
has to be made. This matching has been already discussed in detail in ref. |]13|| . However, 
there we have set q 2 = —m 2 c in evaluating the penguin diagrams of fig. 0. Meanwhile, 
we have realized that in the effective 3-quark theory, q being related to external light 
quarks must be much smaller than m c and the consistent procedure within the operator 
product expansion is to set it identically to zero in the full expression. Numerically, the 
difference between coefficients calculated with q 2 = —m 2 and q 2 = is small, however, the 
resulting modified expressions represent this time correct matching. Proceeding this way 
and applying an analogous procedure to QED penguin diagrams we are able to calculate 
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z(m c ) of ( |4.6| ). We have in the NDR scheme 



^DR ( 



m r 



„NDR, 



y NDR, 



-a s /(24vr)F s NDR (m c 
a s /(87r)F s NDR (m c ) 

-a s /(24vr)F s NDR (m c 
a s /(87r)F s NDR (m c ) 
a/(6vr)F e NDR (m c ) 


a/(6vr)F e NDR (m c ) 




(4.29) 



where 



^NDR 



(A*) 



,m 



ln(-f ) + 1 1 (3z 1 (fi)+z 2 ( f ,)) 



(4.30) 



(4.31) 



with « 0(m c 



In the NDR and HV schemes the "1" in ([OOp and (|OID is absent. 

for z 7^ 1, 2 but in the case 



Consequently, in the NDR and HV schemes one has Zi(m c 
of the NDR scheme these coefficients are non-vanishing. 

Due to the fact that Z\ and z 2 have opposite sign and |zi(m c )| ~ ||z2(?n c )| the coeffi- 
cients Zf(m c ) and Zg(m c ) are very small. This is reminiscent of the situation in exclusive 
two-body decays such as D° — > n°K , for which such a cancellation between Qi and Q 2 



contributions takes place. In the context of the AS 
has already been pointed out in ref. 



1 Hamiltonian such a cancellation 
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Again, when using the formulae above one has to take care that O (a 2 ) terms are omitted 
as well as terms 0(a^) and 0(aa s ) 

At this stage it is mandatory for us to compare the results in (fOO) and ( [4.31 ) with the 
ones of 



30] and where the effects of incomplete GIM mechanism resulting from 



m u ^ m c for /i > m c has been calculated. These analyses resulted in a scheme independent 



constant "|" in place of "1" in ( f4.30| ) and ( ^4.31| ). Yet, the proper matching must involve 



a scheme dependent constant. Our interpretation of this difference is as follows. It is 
intuitively clear that an incomplete GIM cancellation in the Q\ — Q c 2 sector must take 
place for /i > m c . Yet, this effect is at the next-to-leading level at which other scheme 
dependent contributions are present. From our analysis it turns out that in the HV and 



NDR schemes the latter effects cancel exactly the term "|" found in [p9 |, [ pQ] so that 
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Zi(fjj) = for fi > m c . In the NDR scheme however Zi(m c ) ^ 0, although as we will see 



also in this schemes the penguin coefficients are smaller than found in ref. [29 



4.3 The Matching Matrix M 

In order to show that a non-trivial matching matrix M has to be included in the evolution^, 
when going from an /-flavour effective theory to an (/ — l)-flavour effective theory, we write 
an amplitude A for some process as follows 

A = (Q f (m)) T C f (m) = (<?/-i(m)) r C,_i(m) , (4.32) 

where m is the threshold scale between the two effective theories in question. Qj, Cf, 
Qf-i, and Cf-i are the operators and coefficients in these two theories, respectively. 
In analogy to (|2.42| ), one finds from fig. 



(Qf{m)} = (i + ^^f s + ^-6r c ) (5/-i(m)> , (4.33) 

where 

5f i = f < (/)-n(/-l). (4.34) 

The matrices 5f receive only contributions from penguin diagrams which for certain oper- 
ator insertions depend on the number of effective flavours. Inserting ( |4.33|) in (|4.32j ), we 
find 

d/_i(m) = M(m) C f (m) , (4.35) 

where 

M(m) = 1 + ^4^^fJ + ^5fJ . (4.36) 

47T 47T 

A simple consequence of the presence of a non-trivial matching matrix at the next-to- 
leading level is the appearance of small discontinuities in the Wilson coefficient functions, 
and consequently, in the matrix elements of certain penguin operators at the thresholds 
between effective theories. This is not surprising because after all for penguin operators 
Qf Qf-it as explicitly seen in ( |2.1|) . An example of such a discontinuity is given also by 



the initial conditions in (|4.29[) , since for fi > m c one has Zi(ji) — (i ^ 1, 2). The matrices 



M(m) are renormalization scheme independent. They are given in appendix A. 

^This has been pointed out by J.-M. Schwarz and has been briefly discussed in an addendum to E 
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4.4 Master Formulae for Wilson Coefficients (m c < fi < mt>) 

If fi is chosen above m c only operators Qi,2 and Q\ 2 are present in the u — c sector. For 



Qi, i — 1, . . . , 10, we simply find the Hamiltonian in (|4.4j ) with 



^ ] = U 4 (n,m h ,a) M(m b ) ?7 5 (m b ,M w ,a) f Zl |^ W | | ^ ( 1.37) 

v(jj) = U±(p, m b , a) M(m h ) U 5 (m h , M w , a) v(M w ) , (4.38) 
and Zi = for i ^ 1,2. Since 2 enter the effective Hamiltonian with the opposite sign 



to Qi 2 and are only present in the (1 — r) term in (O), we find 

= ~ Z M > = - ^(aO , (4.39) 

and 

y°(n) = zM, = zziji) ■ (4.40) 
On the other hand it is evident from the formulae above that 

yi(M)=3/2(M) = 0, (4.41) 
for arbitrary //. We note that the coefficients y\ and ^ are large. 

4.5 Relations Between Operators 



The operators Qi given in ( |2.1| ) are linearly dependent for / = 4 or equivalently m c < /j < 
Wb- One has the relations 

Q4 = Q3 + Q2 — Qi + Q2 ~ Qi 1 

Q9 = \{ZQx-Qz) + \qI, (4-42) 
Qvo — Q9 + Qa — Q3 ■ 

For / = 3, i.e., /x < m c when the charm quark is integrated out these relations reduce to 

Q4 — Qs + Q2 — Qi , 

Q9 = ~(3Qi-Q 3 ), (4-43) 
Q10 = Q2 + ^(Qi-Qs) ■ 

It should be however emphasized that these relations have been obtained by performing 
Fierz transformations. Consequently, they cannot be used a priori in regularization or 
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renormalization schemes in which the operators Qi and their Fierz conjugates Qi are not 
equivalent in D ^ 4 dimensions. Thus, although they can be used in the HV scheme, 



they receive additional contributions in the NDR and NDR schemes. An inspection of the 
relation (|3.7Q reveals that in the limit a = only the relation for receives additional a s 
corrections so that e.g. in the NDR scheme we find 

Qa = Q 3 + Q2 - Qi - ^ (Qe + Qa - - Iq 5 ) , (4.44) 

which of course then has to be solved for Q4. Yet, we have decided to work with all 
operators and not to use the relations ( f4.42|) -( [4.44|) in evaluating the coefficients C;(//). 



As discussed at length in || one can without any problems work with linearly dependent 
operators. The advantage of such a strategy is a much clearer picture of contributions 
coming from different sources. On the other hand the relations above will turn out to be 
useful in the analysis of hadronic matrix elements. 

4.6 Numerical Results for Wilson Coefficients 

In tabs. 0, [3| we give the numerical values for the coefficients z { and y^ for = 
200, 300, 400 MeV, in the NDR and HV schemes for n = 1 GeV, /i = m c = 1.4 GeV, and 
ji = 2 GeV, at a fixed value of m t = 130 GeV. In tabs. f|, §, [|, we give the m t -dependence 
of the coefficients y-j - yio of the electroweak penguin operators for /i = 1.0 GeV, different 
values of A^jg and the two schemes in question. The coefficients Zi do not depend on 
m t . The m t -dependence of j/3 - y$ is so weak, that it has not been shown explicitly. The 
corresponding results for fi = m c = 1.4 GeV and fi = 2.0 GeV are given in tabs. [7|, ||, |9| 
and 0, [TT], [12] respectively. The reason for giving many tables for several values of /1 is 
as follows: For /i « 1 GeV the calculation of hadronic matrix elements in an approach 
like the 1/N approach can in principle be performed as discussed in the next section. The 
scale /i = 2.0 GeV is very suitable for present lattice calculations. Finally, fi = m c is 
a very convenient scale for the semi-phenomeno logical approach developed in section 6. 
Simultaneously, ji = m c can be used in the case of AC = 1 decays. 

For comparison, we give in the tables mentioned above the corresponding "leading 
order" results. These have been obtained using the strategy of ref. |J. Thus we have set 
all the m t -independent terms in ([4.8|) - (|4.19|) to zercQ and we have used the leading order 



evolution matrix, i.e., one-loop anomalous dimension matrices and one-loop /3-function. 
Clearly, the scale A^o in the leading order expressions cannot be identified with A^g which 

2 )\Ve should warn the reader that the authors of ref. [|l9| keep these scheme dependent terms in their 
leading order evaluation. 
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enters the next-to-leading formulae. In addition, a s at the next-to-leading level is smaller 
than a s in the leading order when A^o — ^ms is taken. In spite of this we have given the 
leading order results for A^o — ^-ms' because the range of values for A^ used in these 
tables overlaps considerably with the range used for the QCD scale in the leading order 
analyses. In any case what is really relevant here are the results for z± and with next-to- 
leading order corrections taken into account, because in this case Aj^g is exactly the scale 



as extracted from next-to-leading order analyses in other processes [J3T 
We make the following observations: 

i) The coefficients z\ and z 2 are suppressed through next-to-leading order corrections 
in both schemes considered, with the effect being stronger in the NDR scheme. Con- 
sequently, the ratio Z-jz+ relevant for the AI = 1/2 rule is smaller than in the 



leading order, in contrast to the statements made in refs. [O, We will clarify 



this separately below. The /^-dependence of z± is shown for Aj^g = 300 MeV in fig. [| 

ii) For /i > m c the coefficients Z{ (i 7^ 1,2) are zero. For /x = m c they vanish in the 
HV scheme and in the NDR scheme they are so small that we do not show them in 
tab. |2|. They remain small in both schemes for /x = 1 GeV although considerable 
enhancements are observed in the case of z^ - Zq for the NDR scheme. The coefficients 
z-j - z w being 0(a), can be neglected for all practical purposes. 

iii) The coefficients y% - y& are very weakly dependent on m t and only their values for 
m t = 130 GeV are shown. The coefficients y§ and y± are larger than y 3 and x/5. We 
notice considerable dependence of y 6 on /x, A^jg and the scheme considered. Whereas 
in the HV scheme y§ is suppressed by ~15% relative to the leading order result, 
in the case of NDR a «10% enhancement of y§ for Aj^g > 300 MeV and /x > m c is 
observed. For /j < m c however, y e in the NDR scheme is slightly suppressed relatively 
to the leading order result although its absolute value is always larger than in the 
HV scheme as shown in fig. ^. The small discontinuity in y 6 for /x = m c in the HV 
scheme is related to the matching discussed in section 4.3. In the NDR scheme this 
discontinuity is larger because in addition z 6 becoming non-zero at /x = m c suppresses 



visibly \y$\ in view of the formula (O). All these changes are compensated by the 
corresponding changes in the hadronic matrix elements and the contributions of other 
operators such as Q\ 2 so that the physical amplitudes remain /x- independent. We 
note that these effects are essentially invisible in the case of y$ as shown in fig. [|. 

iv) The coefficients y-j and y$ show strong m t -dependence as illustrated in fig. |5| and in 
the tables. Also as shown in fig. a sizable mt-dependence is seen in the case of y% 
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and yio- As seen in the tables y-j and y$ show only a weak dependence on \i and A^. 
The corresponding dependences of y$ and y w are much stronger. 

We note that y 9 and y w are substantially larger than y 7 and y 8 and essentially do 
not change when going from NDR to HV. y$ is unaffected by next-to-leading order 
corrections, but yio is considerably suppressed. In ref. ||, where Q$ and Qio have 
been eliminated, the effect of yg and yio has been transferred to yi and y2 which 
vanish in our case. In tabs. || |9] and [12] we give only the m t dependence of the sum 



2/9 + 2/io which is relevant for e'/e. We observe that for m t < 150 GeV the sum 
1 2/9 + 2/io I is enhanced by roughly 20% over its leading order value but for higher 
values of mt this enhancement is smaller. 

vi) The coefficient y-j is essentially the same in the schemes considered. As noticed in 
[|, y-j vanishes for m t ~ 145 GeV in the leading order, and becomes positive for 
larger m t -values. The inclusion of next-to-leading order corrections shifts this zero- 
point to higher values of m t , such as m t m 190 GeV for both schemes. Consequently, 
y 7 remains negative for m t m (9(150 GeV), and is roughly of the same order of 
magnitude as y%. Hence for these values of m% it is substantially enhanced over its 
leading order value. However, for mt ~ (9(200 GeV) the coefficient y 7 is so small 
in both schemes that the contribution of to e' / e can be safely neglected for such 
high values of m t . 

vii) Most interesting is however the enhancement of j/g through next-to- leading corrections 
as noticed already in ref. |19| , where only the HV scheme has been considered. We 



find that this enhancement is smaller in the NDR scheme. We also note that with 
increasing mt the enhancement of yg over its leading order values becomes smaller so 
that for m t ~ (9(200 GeV) the NDR result is very close to its leading order value. 

viii) As a result of different impacts of next-to-leading order corrections on y§ and y§ in the 
schemes considered, there is a strong scheme dependence in the ratio (2/7/3 + y$)/y6, 
which roughly measures the relative importance of the dominant electroweak and 
QCD-penguins in the evaluation of e'/e. We show this in table [13] for = m c . We 



observe a large enhancement of this ratio in the HV scheme which however decreases 
with increasing m t . In the NDR scheme this enhancement is substantially smaller 
and for m t w 190 GeV the LO and NDR results are very close to each other. 



ix) The coefficients in the NDR scheme are either equal to the NDR coefficients or only 
slightly different from them. The largest differences are found for z 6 at fi < m c and 
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for at fi > m c . This has to be kept in mind. 



4.7 Comparison with refs. [12], [13] and [16] 



In refs. p0| , the enhancement of Z-/z + from 1 to roughly 3 for // ~ 1 GeV has been 
found by calculating QCD effects in the leading logarithmic approximation. Subsequently, 
it has been stated in refs. |16| and |12|] that z_/z + is further enhanced by next-to-leading 



logarithmic corrections. Although the calculations in |L2| and [16[ are certainly correct, 
the coefficients on which this conclusion was based, are really not the coefficients z± of the 
operators 0±, but the scheme independent coefficients given by 



z±(n) = z±(n) 



1 - ^B ± 



4tt 



which correspond to 



<o ± 0*)> = <o ± Gu)> 



i 



47T 



(4.45) 



(4.46) 



The coefficients B± for the schemes considered in [|16| (DRED) and ||T2| (NDR,HV), 
are as follows 



(4.47) 



The removal of the 0(a s ) corrections in ( 4.45| ) and (4.46) changes the conclusions 
reached in |L2], |i~6| , so that agreement with our paper is obtained. As stressed already 
in 
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12[ , the procedure of |]12| , |Tq1 can be used, but then also the calculation of hadronic 



matrix elements should include the compensating factors given in (|4.46|) so that in both 
treatments the same amplitude is obtained. 

At this stage it should be recalled that similar correction factors are included in the 
case of the parameter Bk in the QCD analysis of the ^-parameter |3~2]] . However, there 
these corrections are so small that it is irrelevant how they are treated. In the case at hand 
the corrections are larger, and different conclusions about the signs of next-to-leading 
order corrections can be reached, dependently whether z± or z± are considered. Similar 
comments apply to our paper of ref. [Jl3~| , where scheme independent coefficients Z\ - z Q 
and y\ - y 6 have been used. 

Although the introduction of these additional corrections may in principle identify 



the full physical effect, as discussed in [0, it seems to us at present that for the future 
treatment of hadronic matrix elements it is more elegant, more convenient, and more 
transparent not to perform additional "rotations" to scheme independent bases. Indeed, 
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on one hand such a procedure is not unique and on the other hand it may in certain 
cases cause misleading conclusions. The observed suppression of z^/z + through next-to- 
leading order corrections in the scheme considered must of course be compensated by the 
corresponding effects in the hadronic matrix elements, in order to explain the AI = 1/2 
rule as we will discuss in more detail in section 6. 



5 Hadronic Matrix Elements 
5.1 General Remarks 

An important ingredient of any analysis of non-leptonic i^-decays are the hadronic matrix 
elements of operators Qi which we denote by 

(Qi)i = ((tttt)/ \Qi\ K) , 1 = 0,2. (5.1) 

These matrix elements depend generally on the scale \x and on the renormalization scheme 
used for the operators. These two dependences are cancelled by those present in the 
coefficients Ci(fi) so that the effective Hamiltonian and the resulting amplitudes do not 
depend on \i and on the scheme used to renormalize the operators. 

At this point, it should be emphasized that eq. ( |1.2j ) is valid beyond perturbation theory 
so that in principle \i can be chosen completely arbitrary. It can even be set to /i = 0. In 
the spirit of ref. ||33|| , we can then divide the renormalization group evolution into a short 



distance evolution from M w to /i described by 

C(fi) = U(fi, M w , a) C(M W ) , (5.2) 
and a long distance evolution from scale up to fi, described by 

(Q T (fi)) = (Q T (0))U{0,fM,a), (5.3) 

where U(0,/i, a) is the evolution matrix in the long distance regime. Although this evo- 
lution cannot be evaluated at the lower end in a perturbative framework, it has general 
features similar to the evolution in the short distance regime as demonstrated in [R3J. In 
particular, the fi and scheme dependences at the upper end of the evolution (around n) 
must match the ones present in the coefficients Ci(n) in such a way that these dependences 
are not present in the physical amplitudes. In other words it is a matter of choice what 
belongs to the matrix element and what to the coefficient function. 

From this discussion it is evident that the scheme dependence of (Qi{^)) is only present 
at the upper end of the long distance evolution. Consequently, for sufficiently high \i this 
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dependence can be calculated together with the /^-dependence in a perturbative framework. 
The relation between the matrix elements calculated in two different schemes is given in 
fl3.7| ). The relation between matrix elements evaluated at two scales mi < 1712 is simply 
given by 

(Q T (m 2 )) = (Q T ( mi )) U( mi , ma, a) , (5.4) 

where U is the evolution matrix of section 2.3. Thus, if there is a method to calculate 
at a single scale \i — mi, the renormalization group will give us these matrix elements at 
any other scale fi 7^ mi. Unfortunately, the full long distance evolution involves all scales 
down to fi = and consequently, the actual evaluation of (Qi(mi)) can only be done 
in a non-perturbative framework. In spite of this, it turns out to be very instructive to 
analyze the evolution of (Q(/j,)) given in ( |5.4| ) in a range of fi for which U(mi,m2,a) can 
be calculated as in section 4. Since this has not been studied so far in the literature let us 
recall that whereas the evolution of C(fi) in ( |5.2|) in governed by 7 T , the evolution of (Q(/u)) 
is determined by 7. Now the matrix 7 is rather asymmetric ||12||-||15||. In particular the 
elements 7^ = 7^2 for i = 3, . . . , 10 vanish. Consequently, the structure of the evolution in 
( |5.4D is rather different from the evolution in (|5.2| ). Thus, whereas the evolution of Ci i2 (/i) 
is unaffected by the presence of penguin contributions, the evolution of (Qi, 2 (/i)) depends 
on the size of (Qi(fi)), i 7^ 1, 2. Conversely, whereas the evolution of Cj(/i), i 7^ 1, 2 depends 
on the size of Ci )2 (/x), the evolution of the matrix elements of penguin operators (Qi(fi)), 
i 7^ 1, 2 is a sole penguin affair. This different structure brings certain surprises as we will 
see in section 6. 

During the last years there have been many attempts to calculate the matrix elements 
(Qi) by using the lattice approach, 1/iV-approach, QCD sum rules, hadronic sum rules, 
effective QCD action ||34j| , chiral perturbation theory |35] and vacuum insertion. Yet, it 
is fair to say that only in the 1/iV-expansion of refs. |33, |36| and in the vacuum insertion 
approach have all matrix elements (Qi) been calculated. Some subset of the (Qi) is also 



known from the lattice method, sum rules and effective QCD action approach of ref. [34 



At this stage it should be stressed that the 1/N- approach, although bearing some 
similarities to the vacuum insertion (factorization) method, is really a systematic non- 
perturbative method for QCD calculations, whereas this cannot be said about the vacuum 
insertion approach. Moreover, there are important quantitative differences between these 
two methods in the (Qi, Q 2 ) sector. Whereas the 1/N approach offers a plausible descrip- 
tion of CP- conserving K — > tttt amplitudes such a description is not possible in the vacuum 
insertion approach. In particular, the understanding of the AI = 1/2 rule is completely 
missing in the factorization approach. Next, one should mention a complete failure of the 
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factorization method in a large class of exclusive non-leptonic B- and .D-decays, whereas 
a reasonable description of these decays can be achieved in the 1/iV-approach |37, j38| . 
Finally, the value for B^ w 0.70 ± 0.10 obtained in the 1/iV-approach is rather close to 
the present "world average" value given in ( [7.7|) which should be contrasted with B^ = 1 
obtained in the vacuum insertion approach. 

Concerning the matrix elements of the dominant penguin operators the structure of 
the 1/N and the vacuum insertion formulae is similar. However, it should be stressed 
that whereas the scale A x entering these formulae can in the 1/N method be related 
to the ratio F^/F^ |39|, [36], no such relation exists in the vacuum insertion approach. 
Consequently, a quantitative estimate of these matrix elements in the latter approach is 
essentially impossible. 

Concerning the /x-dependence it is well known that the matrix elements of (Qi,Q2) 
calculated in the vacuum insertion method are /x-independent in contradiction with one 
of the general properties of the hadronic matrix elements of operators carrying anomalous 
dimensions. In the 1/iV-approach of ref. the /x-dependence in {Qi(fi)), i — 1,2 arises 
through loop effects in a meson theory and at the semi-quantitative level the matching in \x 
between long and short distance calculations is indeed possible. A problem exists however 



because in the approach of ref. the long distance evolution in ( p.3| ) can only be extended 
to scales /x ~ 0.6 GeV at which perturbative evaluation of the Wilson coefficients C,(/i) 
is questionable. There is a hope that by including vector meson contributions to the long 
distance evolution the scale /x could be increased to /x ~ 1 GeV P0 |. However, such a 
calculation is very difficult and in any case not available for K — > at present. 

On the other hand, the /x-dependence of the for K — > tttt dominant penguin operators 
{Qq{^)) and (Q 8 (/x)) is to a large extend under control in both approaches. It is to a very 
good approximation given by the /x-dependence of the running strange quark mass l/m^/i) 
as we will demonstrate below. This /t-dependence cancels to a very good approximation 
the /t-dependence of C 6 (/x) and Cg(/x) so that in quantities in which Q 6 and Q% dominate, 
the left-over /x-dependence is rather weak. The actual issue then however is what the value 
of m s at a fixed value of /x, say /x ~ 1 GeV is. We will return to this question as well as to 
the /x-dependence of other penguin operators below. 

Next, the renormalization scheme dependence of (Qi(ji)) should be addressed. The 
vacuum insertion method is certainly insensitive to this dependence. It appears at present 
that this is also the case of the 1/N approach. Evidently, in these two approaches the 
scheme dependence of Cj(/x) cannot be cancelled and this feature introduces an additional 
theoretical uncertainty. We will investigate the size of this uncertainty below. 

At this stage a few remarks about the lattice approach are in order [51]- [3B|. In 



29 



this approach the //-dependence and the scheme dependence of (Qi(^)) is in principle 
calculable ||7]]- |52| . However, the accuracy of the lattice methods is insufficient at present 
to contribute very much to the issue of the [i- and scheme dependence at a quantitative 
level. This is evident from the strategy of ref. [19| which we will discuss in more detail 
below. Moreover, the matrix elements (Qi)o and (Q2)o have not been calculated on the 
lattice yet. 

It will be useful in what follows to recall the matrix elements used in refs. || [H], [L9| . 
and to compare them with each other. In section 6, we will go beyond these approaches 
to find matrix elements which we believe are closer to the true QCD matrix elements than 
these given by the methods just discussed. 



5.2 Explicit Formulae for the Matrix Elements 



Here, we will give explicit formulae for (Qi)i which are general enough that the matrix 
elements used in refs. || |7|, [18], [19| could be compared with each other. In order to 
simplify the presentation, we will use the relations between operators given in section 4.5. 
As discussed there care must be taken when using these relations in the NDR and NDR 
schemes. In section 6, we will outline our strategy in using the formulae given below. 
Taking first fi < m c , we have 



(Qi)o 
(Q 2 )o 

(Qs)o 

(Qsh 

(Qe)o 

(Q 7 )o 

(Qs)a 

(Q 9 >Q 

(Qio)o 



1 



-XB\ 
9 1 

XBi 1/2) 



(1/2) 



XB. 



(1/2) 



(Q 3 )o + <Q 2 )o - (<2i)o 
-B^ 1/2 \Q 6 ) , 



ml 



m s (ji) +m d (fi) 



F^ B {i/2) 



g<Qe)o(« + l) - y 



B 



K 



(1/2) 



i: 



X 
~6 



B 



(1/2) 



!(Qi)o-7t(Q 3 )o 



(Q 2 )o + ^(Qi)o-^(Q3)o, 



(5.5 

(5.6 

(5.7 
(5.8 
(5.9 

(5.10 

(5.11 
(5.12 
(5.13 
(5.14 



30 



where 
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X = ^/|^ (m 2 , - ml) = 3.71 ■ 1(T 2 GeV 3 
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For the ease of the reader, we give all the actual numerical values for A x , mg, m n , F n 
that produce these "magic numbers" above in appendix C. 

The overall normalization in eqs. (|5.5|)- fl5T22"D agrees with 0, but differs by a factor 



3/2 relative to the normalization used in |TJ|. This difference cancels however in the final 
result for e'/e and other physical quantities. 

For m c < n < m h the relations ( |4.42| ) are valid and the formulae for (Q^q, (Qg)o and 
(Qio)o change: 



(Qa)q = (Q3)o + (Q2)o-(Qi)o + (Q c 2 )o-(QT)o 

(Q 9 } = \(Qi)*-\(Qz)* + \(q%, 

(Qio)o = (Q 4 )o+ (Qg)o- (Qs)o- 



(5.23) 
(5.24) 
(5.25) 



The matrix elements (Q\ 2 )o,2 although vanishing in the vacuum insertion method may be 
non-zero in general as discussed below. Also eq. (|5.16|) is modified by 0(oi) corrections. In 
the vacuum insertion method B; = 1. 



5.3 Critical Comparison of Various Matrix Elements Estimates 



Having formulae ( |5.5| )-( |5T25|) at hand, we can compare these matrix elements with the ones 



used in refs. 



p, m P 
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In ref. §, eqs. (^5|)-(^22|) have been used with the terms "1" and X in ( PH ), fl5TT2D 
and X in Q5.17Q , ( |5.18|) omitted and B\ AI ^ = 1 except for 



(1/2) 



5.2 



5. 



(1/2) 



2.2 



(3/2) 



0.55 



(5.26) 



in accordance with the 1/iV-approach of ref. [33]. The fact that B^ 2 ' > 1 for i = 1,2 and 
; < 1 allows to come closer to the A J = 1/2 rule than it is possible in the vacuum 
insertion method. Furthermore, fi = 1 GeV and 125 MeV < m s (l GeV) < 200 MeV have 



been used. At this point we should mention that the approach of Pich and de Rafael pi | 
gives similar results to ]33| for K —>■ tttt amplitudes. 

The Dortmund group [|17|, pl| uses essentially the above formulae with 0.6 GeV < fx < 
1 GeV. The main difference with respect to |J is the modified factor B^ 2 ^ 1 for the operator 
Qq. Using the tables of ref. [[17], [18|], we find 



B, 
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2.4 ±0.2. 



(5.27) 



The authors of [[17], [Ij§ also use 100 MeV < m s (l GeV) < 175 MeV. The difference of 



Bq 1 from unity follows according to ref. |17|, [18| from next-to- leading 1/N corrections for 



which the details are unfortunately not available. B$ 
very close to unity. 
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extracted from n7\, 18 



remains 



The Rome group Jl9| uses the set of matrix elements relevant for fi = 2 GeV with 
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(5.28) 



where entries with W are educated guesses and the remaining values are taken from lattice 
calculations [^T|- [1B|. Here, B^ 2 enters the matrix elements (Qi 2 )o as follows 
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(5.29) 



where we have introduced the minus sign in order to make B\ positive. In the vacuum 
insertion approximation B\ 2 = holds. However, in QCD it could certainly be non-zero. 
The range of values for B{ 2 in ( |5.28j ) is an educated guess by the authors of ref. [|19j . As we 
will demonstrate in section 6, the size of {Q\ 2 )o is related to the size of (Qi)o, £ = 3, . . . , 6, 
and it is possible to go beyond educated guesses and calculate B{ 2 . Since as seen in ( [4.40D 



and stressed in ref. |19[ the coefficients y{ 2 are large the operators Q\ 2 could have in 
principle some impact on e'/e. Our calculations in section 6 show that the educated guess 
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of [I9| overestimates B\ 2 by a factor of four, and consequently, Q\ 2 play only a minor role 
for /i < 2 GeV. A detailed comparison with ref. [19| will be given later on. 
Furthermore, the Rome group uses 

m s (2 GeV) = (170 ± 30) MeV . (5.30) 

which corresponds to 

m s (l GeV) = (225 ± 30) MeV , (5.31) 

and is considerably higher than the values used in the other two papers. Consequently, 
the matrix elements (Q 6 ) and (Q 8 ) evaluated at /i = 1 GeV are smaller than in ref. || and 
substantially smaller than in |L7|, |i~8fl . Equivalently, the matrix elements used by the Rome 
group at /i = 2 GeV are close to the matrix elements used in || where = 1 GeV has been 
taken. But since the coefficients C;(/i) of penguin operators decrease with increasing \i the 
penguin contributions to e'/e in ref. |TI| are bound to be smaller than in the remaining 
two papers. We will return to this point below. 

It should also be remarked that whether the terms "1" and X are kept or dropped in 
the matrix elements of Q7 and Qg is essentially immaterial for e'/e. Although (^7)0, ($7)2 
and (Qs)o are modified by roughly 15%, these matrix elements do not play a considerable 
role in e'/e. On the other hand (Qs)2 which is important is modified by only 2%. 

Next, we would like to point out one additional problem with the matrix elements 
given above. They have been calculated in QCD without taking order 0(a) corrections 
into account. It is a simple matter to convince oneself that all the matrix elements listed in 
( ^.5|) -( p r TT^ ) and ( |5.29| ) receive 0(a) contributions because at 0(a) they mix under renor- 
malization with the operators Q7 and Qg which have 0(1) matrix elements. In particular 
one has 

(Q 3 -6> 2 = 0(a) (Ql 2 ) 2 = 0(a) (5.32) 

to be compared with (|5.16|) and {Q\2)2 = used in ref. fll9| . The 0(a) corrections to 
these matrix elements are necessary in order to cancel the scheme dependence present in 
the Wilson coefficients as discussed in detail in section 3. The 0(a) corrections to penguin 
matrix elements should be dropped because they contribute (D(a a s ) or 0(a 2 ) effects to 
the amplitudes. On the other hand the 0(a) corrections to the matrix elements of Qi 
and Q2 should be taken into account because they are exactly at the same level as the 
0(a) effects which have been taken into account when evaluating the initial conditions at 
/i ~ M w in section 4. Similarly O(o) effects to (Qi 2 )o an d (^12)2 have to be included. 

After this rather critical presentation of the existing estimates of hadronic matrix ele- 
ments, we will now turn to an approach which overcomes several difficulties of the methods 
just discussed. 
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6 A Phenomenological Approach to Hadronic Ma- 
trix Elements 



6.1 General Remarks 

We have seen that the present estimates of hadronic matrix elements by means of existing 
non-perturbative methods suffer from many deficiencies, in particular: 

• The hadronic matrix elements of the dominant penguin operators Q$ and C/§ are 
poorly known due to the poor knowledge of ra s (l GeV). 

• The size of hadronic matrix elements of C/i and Q2 is not fully in accordance with 
the known data for K — > ixix . 

• The renormalization and /i-dependences of the hadronic matrix elements are not 
under control. 



These deficiencies imply another weak point of the e'/e analyses of refs. || 17, TE\. The 



usual practice is to vary m s (l GeV) in eq. ( |5.10| ), while keeping the matrix elements 
of current- current operators fixed. Yet, such a procedure fails to fit the data on CP- 
conserving K — > mi amplitudes since for a given value of m s only certain matrix elements 
of current-current operators can fit these data. 

We would like to propose here a phenomenological approach to obtain the hadronic ma- 
trix elements which incorporates properly the renormalization scheme and //-dependences 
of the hadronic matrix elements as given by QCD, and at the same time assures that the 
theoretical CP-conserving K —>■ ixix amplitudes are in agreement with the data. The basic 
ingredients in this approach are 

• the assumption that 7i C Q of ( |4.4|) properly describes the existing data on CP-conservinj 
K — > 7T7T amplitudes. 

• the Wilson coefficient functions calculated reliably in QCD and QED in section 4. 

• some relations between matrix elements which are valid in QCD and are consistent 
with the renormalization group properties. 

Our phenomenological approach to the hadronic matrix elements involves four basic 
parameters: 



the QCD scale parameter in the MS renormalization scheme, 
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• the 5-factors B^ 2 \m c ), B^ 2 \m c ), and B^ 2 \m c ), parameterizing the matrix el- 
ements (Q 2 (m c )) Q , (Qe( m c))o; an d {Q%(j n c))2 respectively, as defined in eqs. (|5~6| ), 
( |5TTD|) , and ( §7Tj ). 

The choice [i = m c turns out to be very convenient as we will see soon, but is not necessary. 

For given values of these four parameters and the chosen renormalization scheme for 
operators our semi-phenomenological approach gives all the hadronic matrix elements of 
the dominant operators while being consistent with the measured CP-conserving K — > ttti 
amplitudes. In this way an interesting correlation between matrix elements of current- 
current and penguin operators absent in the analyses of refs. || [T7|, [Ll| |H| arises. 



Having calculated all matrix elements this way, we will be able to give predictions for 
s'/s as functions of the four basic parameters listed above and of m t . 

In what follows, we will first give the basic formulae of our phenomenological approach. 
Subsequently, we will calculate the hadronic matrix elements (Qi) 0,2 and we will compare 
our results with those discussed in the previous section. 

In our numerical analysis we use the following experimental data 



ReA) = 33.3 ■ 1(T 8 GeV , ReA 2 = 1.50 ■ 1(T 8 GeV , (6.1) 

which express the AI = 1/2 rule 

^-- = 22.2. (6.2) 
6.2 (Qi(fi)) 2 for (V-A)®(V- A) Operators 

It turns out that to a very good approximation all the matrix elements (Qi(A t ))2 can be 
determined from Re^ in ( |6.1| ) as functions of A^jg, fi and the renormalization scheme 
considered. 

We first introduce the Q± operators 

Q± = 7^2 ± Qi) , z± = z 2 ± z\ , (6.3) 

and using ( (4.11) , we find 

ReA 2 = cz + {fi)(Q + {fi)) 2 , (6.4) 

where 

c = ^ V ud V: s = 1.77 • 10- 6 GeV" 2 . (6.5) 

In obtaining (|6.4|), we have set r yi = in view of Zj >> r y h and we have neglected 
the contributions of the electroweak penguin operators which having Z{ ~ 0(a) play only a 
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secondary role in CP-conserving amplitudes. Recall also that {Qi)2 = 0{a) for i = 3, . . . , 6 
and also {Q\ 2 )^ — 0(a). Consequently, ( |6.4j) is valid for any scale /i in the limit a = 0. 
Moreover for \x = m c the contributions of penguin operators and of Q\ 2 (see below) is either 
zero (HV,LO) or very small (NDR) and at this scale (|6.4j) is an excellent approximation 
even for a^O. 

We note next that (see ( f5.15| )) 



(QMh = (g 2 (/x)> 2 = (QMh ■ (6.6) 

This relation is valid in pure QCD in which isospin is conserved. It is consistent with 
the QCD renormalization group evolution. This follows from the fact that the anomalous 
dimension matrix % in the (Qi,Q2) sector is symmetric with equal diagonal entries and 
the fact that {Qi)i = for i = 3, ... ,6. Here, we neglect the mixing with electroweak 
penguin operators which is 0(a). The effect of this mixing will be analyzed below. 

Using the experimental value for Re/^, together with ( |6.4| ) and ( |6.6j ), we find the matrix 
elements ($1,2)2, 

(Ql(/i)>2=(Q 2 (/i)>2 = ^V (6.7) 

Since z + (/i) depends on \i and the renormalization scheme used, ( |6.7|) gives us automatically 
the scheme and /^-dependence of the matrix elements in question. It also gives the Aj^g 
dependence of {Q 1,2)2- 

Using the relation (|5.19|) and ( |6.7|) , we find the matrix elements ((59,10)2 

(QMh = (QiMh = . (6.8) 

As discussed in section 4.5, eq. ( |5.19| ) is valid in all schemes for a = 0, and so is ( |6.8j ) 
in this limit. 

A numerical analysis shows that the /i-dependence in (|6.8[ ) does not cancel the in- 
dependence of Zg(fj,) and zx (fi). This is not surprising, because due to the mixing of Qg and 
Qio with other operators the cancellation of the //-dependence involves more operators. It 
should be stressed that the matrix elements in (|6.7| ) and ( |6.8| ) do not involve any parameters 
except for Aj^g. Although the impact of 0(a) corrections to the determination of these 
matrix elements is very small we would like to discuss here the matrix elements ((^12)2 
which as stated in ( J5.32Q are 0(a). The point is that these matrix elements can be 
calculated ! First they can only be non-zero for /t > m c because for /t < m c the operators 
Ql 2 disappear from the effective theory. Requiring the //-independence of ReA 2 and noting 
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that in the HV and NDR schemes Zi(m c ) = for % = 3, . . . , 10 we find the following "initial 
conditions" 

(Qi(m c )} 2 = (Q c 2 (m c )} 2 = (6.9) 

valid in these schemes. The evolution Q5.4J) of these matrix elements to scales /x > m c 
generates calculable non- vanishing values of 0(a). Using ( |3.7| ) one can also find the cor- 
responding initial conditions in the NDR scheme. The inclusion of these operators in 
necessary if one want to keep KeA 2 /i- independent. These points will be illustrated more 
explicitly in the case of (Qi(/x))o- 

Comparing ( |6.7| ) and ( |6.8| ) with the general formulae ( |5.15| ) and (|5.19|) , respectively, 
we can extract B^ 2 \fj) which with the approximation used is common for Qi, Q 2 , Q 9 
and Qi . The result is represented by the solid lines in figs. [FJ and |8] for Q 12 and Q^xo, 
respectively. The remaining lines in these figures represent exact RG evolution. They are 
obtained by evaluating first the matrix elements at /x = m c with the help of (|6.7|) and 
and evolving subsequently to other values of /x by means of the evolution equation 
In the evolution we have used A^g = 300 MeV. 

We observe the following points: 



Our formulae ( 577 ) and ( |5.8| ) are good approximations to the exact evolution including 



0(a) effects for 1 GeV < /x < 2 GeV, considered mainly in this paper. For (Qi) 2 
and {Qg} 2 the approximations become slightly worse when /x is increased beyond this 
range. 

We observe a sizable /x-dependence of b[ 3 ^ 2 \/j,). 

The next-to-leading order corrections to Zx$(fi) decrease B^ (/x) relative to the 
leading order determination with the effect being stronger in the case of the NDR 
scheme. 

Most importantly however, we find B^y 2 \/j, fa 2 GeV) fa 0.44, to be compared with 
B^ 2 ^ fa 0.8 ± 0.2 in ( 5.28 ). We conclude therefore that the lattice estimate for 



(Q9,io(2 GeV)) of ref. [|1^] overestimates these matrix elements by almost a factor of 
two. We will investigate the impact of this on e'/e in section 8. The 1/N approach 
used in is in a much better shape here because for low /x, b[ 3 ^ 2 ^ fa 0.55 in ( |5.26| ) 
is quite consistent with figs. [7J and |8| for /i fa 1 GeV. 



37 



6.3 (Qi(/i)) for (V - A) <g> (V - A) Operators 

The determination of {Qi{jj))o operators is more involved because several operators may 



contribute to ReA and a relation like Q6.6Q does not exist in this case. Yet, some progress 
beyond the approaches of section 5 can also be made in this case. As we will see the choice 
H = m c is again very helpful in this respect. 

Let us begin with \i > m c . In this case the Wilson coefficients Zj of penguin operators 
are zero and ReA$ is given as follows 



ReAr 



z-(jJ.) (<Q-)o - (Q-)o) + z + (ji) ((Q + ) Q - (Ql)o)} , (6.10) 



where similarly to (|6.4|) , r jji have been set to zero. On the other hand, for fi < m c the 
operators Q± disappear from the effective theory but instead the QCD penguin operators 
contribute 

ReA = c [z_(im) (Q-(j*)) + z+(jj,)(Q+(im)) ] + Re< , (6.11) 

where 

6 

Re< = c^2zi(n) (Qi(jz))o. (6.12) 

i=3 

We again neglect the contributions of electroweak penguin operators Qi, i — 7, . . . , 10. 
Since ReA Q cannot depend on /i, expressions (|6.10j) and (]6.11|) have to be matched 



properly at /i = m c , the transition point at which one goes from / = 4 to / = 3 effective 
theory. This matching gives the following relation between (Q±(m c ))o and (Qi(m c ))o, 
i = 3, . . . , 6 operators 

6 

^Zi(m c ) (Qi{m c )) = - z-{m c ) {Q c _(m c )) + z + (m c ) (Q+{m c )) , (6.13) 

i=3 

which clearly shows that the role of penguin operators in the AI = 1/2 rule is for \i > m c 
played by the operators Q\ 2 . There is a very immediate consequence of this observation. 
The size of the matrix elements (Q±)o is correlated with the size of (Qi)o, i = 3, . . . , 6. 

Now, as discussed in section 4, in the HV and NDR schemes Zi(m c ) = holds for 
i = 3, . . . ,6. Consequently, in these schemes we have the relation 

2-W (Q-(m c )) + z + {m c ) (Q c + (m c )) = , (6.14) 

which is also valid for 

If we next make a very plausible assumption consistent with all existing non-perturbative 
methods that (Q±) > 0, we find the following "initial" conditions for the matrix elements 

of Qy 

(QiK))o= (Q c 2 K))o = 0, (6.15) 
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valid in the HV and NDR schemes. Using ( |3.7| ) with a = 0, we then find the corresponding 
result for the NDR scheme 



(Qt(m c )) = 0, (6.16) 
aJm c ) 



(Q c 2 (rn c )) 



12vr 



(Q 6 (m c )) + (Qi(m c )) - ^(Q 3 (m c )) - ^(Q 5 (m c )) 



One can easily check that this result is consistent with Q6.13 ) and with the initial conditions 
for Zi(m c ), i = 3, . . . , 6, given for the NDR scheme in ( [4.29 ). The effect of 0(a) corrections 
can easily be taken into account by using (ft.7|) . It has been included in our numerical 
analysis. With the initial conditions ( |6.15|) and (|6.16|) at hand we can calculate (Qi(^))o 
for any /i > m c . The result is given in section 6.5. 

We now want to constrain the matrix elements by the experimental value of ReAo. To 
this end it is useful to set \x = m c in (|6.11|) . Since in the HV scheme for /i = m c only Qi 
and Q 2 operators contribute to ReAo, we find the relation 

z\{m c ) (Qi(m c )) + z 2 (m c ) (Q 2 (m c )) = , (6-17) 

from which the following expression for the matrix element (Qi{m c ))o as a function of 
(Q 2 (m c )) Q can be found 

(Qi (m c = - f r : ^{Q 2 (m c ))o, (6.18) 

Using next the relations ( |5.8| ), ( |5.13| ), and ( |5.14| ), we are able to obtain (Q4,(m c )) , (Q9(m c )) , 
and (Qio(^c))o- Because (Q3(m c ))o is colour suppressed, we could set it to zero or just 
use the expression of eq. (|5.7|). Since Q3 has small Wilson coefficients it does not play any 
role in our analysis. On the other hand (Qi{m c ))Q plays a substantial role in the anal- 
ysis of e'/e and the assumptions about (Q3(m c ))o entering (|5.8|) matter to some extent. 
In the following, we will therefore take (Q3(m c ))o according to eq. (|5.7|). (<5g(?77, c )) and 
(<3io( m c))o are much less important for e'/e than (Q4(m c ))o. 

In the case of the NDR scheme there is a small contribution of the penguin operators 
to ( |6.15| ). Moreover, the determination of (Qi{m c ))o requires the use of the relation ( |4.44| ). 
Thus, we need some assumptions about (V — A) ® (V + A) operators which are discussed 
below. For \x = m c however the issue of (V — A) <g) (V + A) operators has only very 
small impact on fl6.18|) even in the NDR scheme. Equivalently, the matrix elements of 
(V — A) (gi (V — A) operators in the NDR scheme can be obtained from the HV matrix 
elements by means of 
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With all this formalism at hand, we can now calculate the parameters (m c ) and 
B^' 2 \m c ) as functions of B2 (m c ). If we in addition make the very plausible assump- 
tion valid in all known non-perturbative approaches that (Q-(m c )) > (Q + (m c )) > the 
experimental value of ReAo implies 

B$3(mc) = 5.8±l.l, B^ K (m c ) = 6.7 ± 0.9 , B^{m c ) = 6.3 ± 1.0 . 

(6.19) 

This should be compared with the vacuum insertion result B^ = 1, and B^ ~ 2.2 in 
the 1/iV approach. The latter value corresponds to fi ~ 0.6 GeV, and when extrapolated 
to jj, = m c would give -83 (to c ) ~ 2.8. Yet, it is clear that the imposition of the AI =1/2 
rule enhances the role of (V — A) <g> (V — A) operators. Using ( |6.18| ), we plot Bf' 2 '(m c ) 
in fig. I as a function of B^ 2 \m c ). We note that if the AI = 1/2 rule is imposed, the 
parameters B^ 2 \m c ) and B^ 2 ^ (m c ) are strongly correlated. A similar plot for B^ 2 \m c ) 
is given in fig. |TD|. Here, the dependence on B^ 2 \m c ) has been shown. These values are 
consistent with the upper end of the range for B^ 2 ^ used by the Rome group |fL9"| . The 
strong deviation of B^ 2 ^ from unity will have interesting consequences for e'/ e as already 
pointed out in 0, and even stronger emphasized in []nj. 



6.4 (Qi(/i)) ,2 for (V - A) <g> (V + A) Operators 

We have seen that by setting /i = m c , we have decoupled the penguin operators from 
the question of the AI = 1/2 rule. The matrix elements of (V — A) <g> (V — A) penguin 
operators could still be determined due to the operator relations of section 4.5. On the 
other hand the matrix elements of (V — A) <g> (V + A) operators cannot be constrained by 
CP-conserving data unless new relations between (V — A) <g> (V — A) and (V — A) <g> (V + A) 
operators are found. Some relations of this type have been suggested in ref. |L7], [18|, p4[ . 
However, in our opinion these relations are suspect as we will see below and therefore, we 
will not use them here. 

For (Qe( m c))o, we simply use the formula ( |5.10| ), and for {Q5(m c )) the relation (|5.9|) , 
however setting B^ 2 ^ = B^ 2 \ The assumption about (Q5{m c ))o is of little importance 
because this matrix element, similarly to (Q3(m c ))o, is colour suppressed and has very 
small Wilson coefficients. Some support for this strategy will be given soon. 

For the matrix elements of the (V — A) <g> (V+A) electroweak penguin operators Q7 and 
Qs, we will simply use the relations (|5.11|) , (|5.12|) , as well as ( p.!7|) and (|5.18|) , modified by 
setting 5{ 1/2) = Bi 1/2) and Bf /2) = Bf /2) . Since mainly (Q 8 (m c )) 2 is of importance for 
e'/e, the simplification of identifying these 5-parameters is immaterial. Again our analysis 
presented below supports this strategy. 
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At this point, a remark on the strange quark mass is in order. Throughout our analysis, 
we have used the value m s (m c ) = 150 MeV, which corresponds to m s (l GeV) = 175 MeV. 
This value is somewhat lower than the central value cited by the Particle Data Group 



However, the values there have been obtained with lower values of Aj^g than are used in 
this work, and since m s (l GeV) decreases for larger values of A^g, we expect the above 
value of m s to be in the right range []55|]. Nevertheless, due to the existing uncertainty in 
the value of m s , we have allowed for a larger variation of B^ 2 ^ and B^ 2 \ which in this 
way include a possible variation in the strange quark mass. 

As seen in ( |5.28|) the parameters B^ 1 ^ and -B^g calculated in the lattice approach |^1 



46| , agree very well with the vacuum insertion value and in the case of B^ 2 ^ and B^ 2 ' 



with the 1/N approach of ref. |25|, PSR . The question arises however, whether this result is 
valid for a large range of /i. It has been stressed in ref. ||, that the /^-dependence of (<5e(A i ))o 
and {Q%{\x))i is governed in the 1/N approach by the /i-dependence of m s . However, to our 
knowledge it has never been checked whether B^ 2 ^ and B^ 2 ^ are indeed /^-independent 
when /i is varied say in the range 1 GeV < /i < 4 GeV. With the evolution equation 
( |5.4|) at hand, we can answer this question. Taking B^ 2 \m c ) = 1 and B^ 2 \m c ) = 1, we 
have calculated the //-dependence of these parameters in the case of LO, NDR and HV. As 
seen m figs. [TT] and 0, B^ /2) and B$ l2) are only weakly dependent on /t, demonstrating 
that at least for these matrix elements the results of present non-perturbative approaches 
are consistent with the renormalization group evolution. On the other hand, one would 
naively expect that the /i-dependence of the colour suppressed matrix elements (Q s(/i))o 
and {Qi{ij))2 should be quite different from the other two matrix elements implying strong 
/i-dependence of B^ 2 \ After all, the Wilson coefficients y 5 and y 7 have a /i-dependence 
very different from y 6 and y$. This naive expectation turns out to be completely wrong ! 
As shown in figs. [13] and [14] also B^ 2 ^ (/t) and B^ 2 ^ (/t) are only weakly dependent on /i. 
Moreover the ratios B^ 2 ^ (/t) / Bz*~ (/i) and B^ 2 ^ (/t) / B^ (/t) are almost independent of 
/i. In our opinion this is a non-trivial result which gives some support to the lattice results 



of Mi- M 



In this connection it is interesting to investigate whether the relations between various 
operators proposed in refs. fT7|, |TB|j are consistent with the renormalization group evolution. 



Taking the relations of Wu [54| at face value and using the general expressions ( |5.5| )- (|5.19| ), 
one can derive relations between the ^-parameters. The most interesting among these 
relations are 

B [ 1/2) = h - Bi 1/2) , M 1/2) = Bi 1/2) = jj Bi 1/2) , Bf /2) = | B? /2) . (6.20) 
It is easy to verify by using Q5.4p that the relations ( |6.20| ) are incompatible with the 
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renormalization group evolution and can only be valid at a single value of //. The question 
then is at which one ? 



6.5 The Fate of the Operators Q\ s 



The authors of ref. |1| presented an educated guess for the size of the matrix elements 



(<5i 2 (2 GeV))o, and emphasized the possible importance of these operators. Having the 
initial conditions for these matrix elements at // = m c given in ( 6.15| ) and ( |6.16| ) for the 



HV and NDR scheme, respectively, we can verify whether this is indeed true. In fig. [15 
we show Bf 2 (fi) defined in Q5.29Q as functions of B^ 2 \m c ) for ji = 2, 3, and 4 GeV in 



the HV scheme. For \i = 2 GeV and B^ /2 \m c ) = 1.0 ± 0.2, as used in [O, B c 12 turn 



out to be at the lower end of the range in ( 5.28 ). We conclude therefore that at least at 



fi = 2 GeV the role of the operators Q\ 2 is rather minor in spite of large coefficients y\ 2 
given in (}4.40|) . They are more important for /i = 4 GeV, but since for < m c they do 



not contribute at all, the fate and the role of these operators depends rather strongly on 
the choice of fi and also as seen in fig. [15] on B^ 2 ^ . Similar conclusions are reached in the 
NDR scheme. 

6.6 Strategy of the Present Paper 

Making plausible approximations, using the experimental data for CP-conserving ampli- 
tudes, and having calculated the Wilson coefficients in section 4, we were in a position 
to calculate the matrix elements (Qi(/i))o,2 and (Qi.2(/ i ))o,2 for any jj, > 1 GeV, and any 
renormalization scheme in terms of the four basic parameters 

Ams, P4 1/2) K), 4 V2) K), Bf l2 \m c ). (6.21) 

The remaining two S-parameters B^ 2 ^ and -Bg 1 ^, play only a minor role in the analysis 
of e'/e and will be set to 1. 

Let us summarize the basic ingredients of this approach: Having fixed the basic param- 
eters in ( |6.21|) the calculation of the hadronic matrix elements in a given scheme proceeds 



in several steps as follows. 
Step 1: 

(Qi(/-0)2, (^(aOK (Q 9 (a0)2 and (Qio(a0)2 are given directly by ( ^7|) and (|6Tg 



Step 2: 

The initial conditions for matrix elements (Qi(/i))o at fi = m c can be calculated as 
described in section 6.3 with the coefficients Zi(m c ) calculated in section 4. 
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Step 3: 

Matrix elements (Qj(//))o for |i / m c can be found by using the evolution for the 
operators given in ( |5.4j ). This can also be done for (Qi(fi))2 starting at ji — m c where 
and (16.81) are most accurate. 



Step 4: 

The matrix elements (Qi(fi))o,2 in any other scheme can be calculated by using the 
relation (|3~7|). 

In performing this program, care must be taken as in the case of the Wilson coefficients, 
that higher order terms in a s and a generated in steps 3 and 4 are discarded. When this 
is done, it is an easy matter to convince oneself that the resulting physical quantities such 
as A , A 2 and e'/e are independent of \i and the scheme considered. 

Now, whereas B 2 (m c ) and the related {V — A)®{V — A) parameters taken in a given 
scheme can be extracted from CP conserving data as seen in ( |6.19| ) and in figs. |T|^T0| this 



is not the case for B^ 2 \m c ) and Bg 2 \m c ). From ( |3.7|) , we know that if B^ 2 \m c ) = 
Bf ,2 \m c ) = 1 in the HV scheme then for = 300 MeV, we find B ( ^ /2 \m c ) = 0.84 
and B^ 2 \m c ) = 0.85 in the NDR scheme. Yet, since we do not know at present which 
renormalization scheme exactly gives the vacuum insertion results, we decided to proceed in 
the phenomenology of e'/e as follows. For all the cases considered (LO, NDR, HV), we will 
take universal values for the parameters B$ 2 \m c ) and B^ 2 \m c ). This will necessarily 
introduce a scheme dependence in our results for e'/e because the Wilson coefficients are 
scheme dependent. This left-over scheme dependence in e'/e will give us some idea about 
the uncertainty in this ratio resulting from the poor knowledge of the scheme dependence 
for (Qe(m c ))Q and (<5s( m c))2- On the other hand the scheme dependence of the matrix 
elements of the (V — A) (g) (V — A) operators will be properly taken into account by means 
of our approach. 



7 CKM— Parameters and e 

In the analysis of e'/ e we will need the values of various CKM-parameters. We use the stan- 
dard parameterization of the quark mixing matrix |k| in which the four basic parameters 
are 

Sl2 = | Ks | , Sl3=|Kb|, S 2 3 = |Kb| , 8, (7.1) 

with 5 denoting the sole complex phase of this matrix. 
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In our numerical analysis we will set 



WJ = 0.221 



\v uh /v c 



Wc 



cb 



0.043 ± 0.004 



cb 



0.10 ±0.03 



(7.2) 
(7.3) 



and as usual we will extract 5 by fitting the theoretical expression for the parameter e of 
the indirect CP-violation to the data. To this end, we use 



where 



exp(m/4) 
V2AM 

Iim4 



(ImM 12 ±2£ReM 12 ) 



(7.4) 



ReA 

with the amplitude A = A{K — 
diagrams 



AM = 3.5 • 10~ 15 GeV , (7.5) 
7T7r)i = o) and M 12 obtained from the standard box 



M 



G\ 



.1.2 



12 7T 2 



FlB K m K M^ \f Vl S(x c ) + \f V2 S(x t ) + 2X* c X* t r] 3 S(x c ,x t ) 



(7.6) 



Here S(xi) and S(x c ,Xt) are the Inami-Lim functions for which explicit formulae are given 
e.g. in ref. 0. Next, F K = 161 MeV, A; = V^V^ and B K is the renormalization group 
invariant parameter describing the size of the matrix element (K°\(sd)v-A{sd)v-A\K°) . 
In our numerical analysis, we will choose 



B K = 0.65 ±0.15 



(7.7) 



which is in the ball park of most recent estimates |5q| . For the QCD-factors r/i, we will use 



0.85. 



r/2 = 0.58, 



0.36. 



{7.1 



where r/ 2 includes next-to-leading order corrections calculated in ref. P2]| . Because the 
second term in eq. ( |7.6| ) is the dominant one, we included next-to-leading order corrections 
in i]2 although for the other two smaller terms they are still unknown. Without these 
corrections one would have (r/ 2 )i Q w 0.62. 

For fixed values of the parameters in (|?T2|), ( |7.3j ), ( |7.?| ) and (|7.8| ), and for a fixed value 
of m t comparison of (|7.4| ) with the experimental data for e gives two solutions for the 
phase 5 in the first and second quadrant. With the experimental constraint from B° — B°- 
mixing it is sometimes possible to exclude one of the two solutions. In particular for 
m t > 150 GeV and > 200 MeV the solution in the first quadrant is favoured. Due to 
remaining substantial uncertainties in Fb and m t , we will however not use the constraint 
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coming from B° — 5°-mixing and we will present the results for e' / e corresponding to two 
solutions for 5. 

In the left half of fig. [IB], we show ImA t as a function of m t when |Kb|, |Kib/Kb| and 
Bk are varied in the ranges given by ( |7.2| ), ( |7.3| ), and ( |7.7|) . In the right half of fig. [L6j 
similar plots are given for the parameter ranges to be expected in a few years time 

| Kb | = 0.043 ±0.002, \V uh /V ch \ = 0.10 ±0.01, B K = 0.65 ±0.05, (7.9) 

We observe that ImAt decreases with m t and is larger for the solution < 5 < n/2. The 
uncertainty in ImA t is considerable at present. 



8 e'/e Beyond Leading Logarithms 
8.1 Basic Formulae 

The parameter e' is given in terms of the amplitudes Aq and A 2 as follows 

to 



where 



ImA 



-_£(l-fi)exp^, 



ReA 2 _ 1 lmA 2 
to = — — — , \l = — - — — 
ReA us Im^o 



ReA ' 

and <p = 7r /2 + S 2 — S ~ 7r/4. 

Using H c s(AS = 1) and the experimental data for u, ReA and e, we find 

ImAt 



- = 10" 

£ 



1.7- 10" 4 



p(l/2) _ p(3/2) 



where 



p(l/2) 
p(3/2) 



E p(l/ 2 ) 
E p(3/2) 



rEyi(ft)o(i-n 

a; 



77+77' y ; 



with 



1.7 



594 GeV" 



U) 
5.2) 



(8.3) 

(8.4) 
(8.5) 

2 |e I ReA ^'^ 
In (|8.5|) and (|8.5|) the sum runs over all contributing operators. This means for \i > m c also 
the contribution from operators Ql 2 to p( 1//2 ) and p( 3 / 2 ) have to be taken into account. 
This is necessary for pt 1 / 2 ) and p( 3 / 2 ) to be independent of /1. Next 

1 (ImA 2 ) LB . 

u ImAn ' (8 - 7) 







77+77' 
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represents the contribution of the isospin breaking in the quark masses (m u 7^ m^). For 
fi,j+r/, we will take 

Slrj+tf = 0.25 ± 0.05 (8.8) 

which is in the ball park of the values obtained in the 1/N approach |25] and in chiral 
perturbation theory |j7], |2"T|| . fi^+j/ is independent of m t . 



Now the formulation of the phenomenology of e'/e found in refs. and [17], 18, [19 

uses instead of P^ 1 / 2 ) and p( 3 / 2 ) defined in fO) and dOI) the ratios 



1 2 P (1/2) (3/2) P {3/2) 

n i = ~7j(T72y ' Q i = TjJN ( l ~*W)» ( 8 - 9 ) 

with fi^+j/ being treated separately. 

We would like to point out one drawback of this formulation: The separate contributions 

(1/2) (3/2) 

and fir depend on /x and on the renormalization scheme. In this respect the 



example of the operator Q\ is very instructive. As pointed out in ref. ||19|| , Q\ can have 



some impact on e'/e for \i « 0{2 GeV). Yet, for /i < m c the operator Q c 2 is absent in 
the effective theory, hence fi^ 2 ' 1 = 0, an d its effects are hidden in the contributions of 
operators present in the effective three quark theory. 

On the other hand, P^ 1 / 2 ) and p( 3 / 2 ) in (|8.4j ) and (|8.5|) are independent of /i and the 
renormalization scheme considered. We have verified numerically that this is always the 
case. Consequently, they are more suitable for the phenomenological analysis. P^ 1 / 2 ) is 
dominated by the contribution of the Q§ operator but receives also important contributions 
from Q4. p( 3 / 2 ) receives only contributions from electroweak operators Qi (i — 7. . . 10). 
Although Q s dominates p( 3 / 2 ), the operators Q7, Q 9 and Q w have also some impact on 
the final result for e'/e 

8.2 ^-Expansions and the Four Dominant Contributions to e'/e 

The contributions P^ 1 / 2 ) and p( 3 / 2 ) can be written as linear combinations of the Pj- 
parameters introduced in section 5. In our approach of section 6, there are only three 
relevant Pj-parameters introduced in ( |6.21| ). We then find 



P^ = a£ m + 4 V2) 2# /2) + 4 1/2) 4 1/2) , (8.10) 

= 4 3/2) + 4 3/2) Pf /2) , (8.11) 

where the parameters Pj will be taken at /i = m c . Then the coefficients af~ and a\ 3 ^ 2 ^ 
depend only on Aj^g, m t , and the renormalization scheme considered. These dependencies 
are given in tabs. [T3 and ITS . 
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Looking at the expression (|8.1Q) and ( |8.11| ) , we identify four contributions which govern 
the ratio e'/e: 

i) The contribution of (V — A) £g> {V — A) operators to P^ 1 / 2 ) is represented by the 
first two terms in ( |8.10| ). These two terms include also small contributions from 
(V — A) (g) (V + A) electroweak penguin operators, and if we would use \i > m c 
also operators Q\ 2 would contribute here. These two terms are dominated by the 
contribution of the operator C/ 4 . We observe that the sum of the two terms in 
question is negative and only weakly dependent on B% 2 \ especially if we vary this 
parameter in the ranges given in ( |6.17D . The dependence of A^jg is also weak. These 
weak dependences on renormalization scheme, B^ 2 ^ and Aj^g result from the fact 
that in our approach of section 6 the matrix elements entering the first two terms in 
p(i/2) are more or i ess fixed by the experimental value of Aq. The weak dependence 
on m t results from the contributions of electroweak penguin operators. It should be 
stressed that the constraint from CP-conserving data suppresses considerably e'/e. 
In order to see this we calculated the sum of the first two terms in (|8.10| ) using the 
matrix elements with Pj = 1. Taking A^g = 300 MeV, fi = m c and m t = 130 GeV 
we have found a[, 1/2) + a^ 1/2) P^ 1/2) equal to -0.75, -0.72, -0.74 for LO, NDR and HV, 
respectively. This should be compared with -3.78, -3.99 and -3.92 used here. 

ii) The contribution of (V — A) <g> (V + A) QCD penguin operators to 

pd/2) 

is given by 

the last term in ( |S.10|) . This contribution is large and positive. The coefficient ag 1 ^ 
depends sensitively on A^g which is the result of the strong dependence of y G on the 
QCD scale, a^ 2 ^ is suppressed by next-to-leading corrections at /i = m c although 
due to matching effects discussed in section 4 it is somewhat enhanced for ji > m c in 
the NDR scheme. 

iii) The contribution of the (V — A) <g> (V — A) electroweak penguin operators Qg and Qio 
to p( 3 / 2 ) is represented by the first term in p( 3 / 2 ) . As in the case of the contribution 
i), the matrix elements contributing to fixed by the CP-conserving data. 
This time by the amplitude A 2 . Consequently, the scheme dependence and the Aj^g 
dependence of a^ 2 ^ are weak. The sizable m t -dependence of a^ 2 ^ results from 

/2 s ) 

the corresponding dependence of yg + y±o. a contributes positively to e'/e. As 
seen in tab. [15] this contribution is slightly enhanced through next-to-leading order 
corrections. On the other hand it should be stressed that it is considerably suppressed 
by the constraint from CP-conserving data. Indeed as seen in fig.lP^ 2 ) (m c ) is by 
a factor of 2 below the vacuum insertion value. 
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iv) The contribution of the (V — A) <g) (V + A) electroweak penguin operators Qj and 
Qs to p( 3 / 2 ) is represented by the second term in ( |S.ll ). This contribution depends 



sensitively on m t and A^g. It contributes negatively to e'/ e. The m t dependence of 
this contribution is governed roughly by 2/7/3 + ys- For m t < 150 GeV enhancement 
through next-to-leading corrections is observed but for m t > 150 GeV the leading 
and next-to-leading order results are quite similar to each other. For m t < 150 GeV 
the operator Q-j is important when next-to-leading corrections are taken into account 
but for m t « 200 GeV its contribution can be neglected. 

The competition between these four contributions depends on m t and on the values of B^ 2 ^ 
and Pg . Setting first B^ 2 ^ = B$ = 1 , we observe that the contribution ii) is always 
largest and the contribution iii) smallest. For lower values of mt the contribution i) is more 
important than iv), and it is this contribution which represents the main suppression of 
e'/e at m t < 150 GeV together with Q v+V > 7^ 0. With increasing m t , the contribution iv) 
grows and for m t > 170 GeV, it is more important than the contribution i). At such high 
values of m t , the collaboration of i) and iv) is so strong that it becomes as important as 
the contribution ii) and hence, e'/e becomes very small. For even higher values of m t the 
sum i)+iv) wins the competition, and e'/e becomes negative. 

8.3 The ^-Dependence of Various Contributions 



It is instructive to study the //-dependence of the various contributions in fl8.4|) and flOf) . 
In fig. O we show P^ 2 \ P^ 2 \ P^ 2 \ Pg 3//2 ' ) and the remaining contributions to P^ 1 / 2 ) and 



as functions of \i for Aj^ = 300 MeV, m t = 130 GeV and B^ /2 \m c ) = P ( {1/2) ' 



V MS 

B?'(m c ) = Bg 2 \m c ) = 1. We observe 



i) The various contributions show visible /x-dependence. The sum of all contributions 
is however /^-independent. The discontinuities in NDR are the result of the matching 
at m c . 

ii) The /i-dependence of Pg (//) and P 8 (/i) demonstrates the importance of 1/N 
corrections because in the large- N limit they should be //-independent. These /1- 
dependence although visible is rather weak. 

iii) Substantially stronger /i dependence is observed for the remaining smaller contribu- 
tions in accordance with our previous discussions. P^ 2 ' remains however always 
very small. 
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8.4 Final Numerical Results for e'/e 

We are now in the position to present our final results for e'/e which are based on the 
Wilson coefficients of section 4, the hadronic matrix elements of section 6 and results on 
ImA t of section 7. 

Let us first emphasize that with the help of the CP-conserving data we were able to 
determine two out of four contributions to e'/e discussed in section 8.2. These are the con- 
tributions i) and iii). Taken together these two contributions give a negative value for e'/e 
which is very weakly dependent on A^g and on the renormalization scheme (LO,NDR,HV). 
For instance for = 300 MeV and m t = 130 GeV we find e'/e « -2.7 x 10~ 4 . We 
observe that with increasing m t this part of e'/e becomes smaller. 

We should also stress that these two contributions are very different in the vacuum 
insertion method. In fact in the latter method the sum of the contributions i) and iii) gives 
a positive e'/e ~ 1.8 x 10~ 4 for m t = 130 GeV instead of the negative contribution found 
here. 

The fate of e' / e in the Standard Model depends then mainly on the following quantities: 
ImA t , mt, A m , B^ 1/2) (m c ), 4 3/2) (m c ) . (8.12) 



This dependence is shown in tabs. |16] pl| and in figs. |i8|^20|. In obtaining these results 
we have set m s (l GeV) = 175 MeV but we varied the parameters B^ 2 \m c ) and B^ 2 \m c ) 
in a broad range which effectively could take into account variations in m s . Recall that 
the relevant combinations are B^/m 2 . B e8 = 0.75,1,1.25,1.5 and 2.0 corresponds then 
effectively to m s (l GeV) = 202, 175, 157, 143 and 124 MeV if the matrix elements (Q 6>8 ) 



were kept at the vacuum insertion values. The tabs. |16H18| in which Bq = B$ show then 



effectively the variation of m s (l GeV) in the range 143 MeV < m s (l GeV) < 202 MeV. 



Tabs. [K], ^0| and [21] on the other hand show examples of Bq ^ Bg. We observe the following 



main dependencies: 

• e'/e decreases with increasing m t and B. 



(3/2) 
8 



• e' / e increases with increasing A^jg and B^ 2 ^ 
Moreover as one can deduce from the analysis of ImA t 

• e'/e increases with decreasing £>k, Kb and V^b- 

The main conclusion which we can draw on the basis of these tables are as follows 
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i) If Bq = B$ = 1 the values of e'/e are below 10 -3 and fully consistent with the result of 
the experiment E731. The next-to- leading corrections calculated in the NDR scheme 
suppress e'/e by roughly 10-20% relative to the leading order result. Much stronger 
suppression is observed in the HV scheme, especially at higher values of m t . Yet in 
the full range of the parameters considered, e' / e remains positive. 

ii) We stress however that the increase of (Qe(m c ))Q by only a factor of two with 
{Qs{ m c))2 kept fixed moves e'/e as shown in tab. [21] in the ball park of the NA31 
result if m t < 150 GeV, > 300 MeV and the NDR scheme is considered. The 
values in the HV scheme are lower but still consistent with the NA31 result if the 
CKM phase 5 is in the first quadrant. 



iii) For the intermediate choices of the Bi parameters shown in tabs. [17], [H] and |2(] 
the ratio e'/e is found somewhere in between the values of the two experiments in 
question although for higher values of m t and lower values of A^jg they certainly 
favour the E731 experiment. 

iv) An important point to be stressed here is the relatively strong dependence of e'/e on 
the value of A^jg. It is a result of large anomalous dimensions of penguin operators. 



This dependence is evident in the tables but is even more clearly depicted in figs. [18 
19| and ^0] where e' / e is plotted as a function of A^g for various choices of m t and 
Bi. We note that for the cases (B$ = Bg = 1.5) and (B& = 1.5, Bg = 1) with 
m t = 130 GeV the increase of A^ from 150 MeV to 450 MeV moves essentially e'/e 
from the E731 value to the NA31 result. 

We should stress that although our different contributions to e'/e have been presented 
for fi = m c we have checked that the final result for e'/e does not depend on the actual 
choice of \i. 

In summary, our results for e'/e are fully compatible with both experiments when 
m t < 150 GeV but favour E731 when m t is higher. The progress in predicting a more 
accurate value of e'/e can only be achieved through the reduction of the uncertainties in 
the values of the parameters listed in ( fg.l2j ). 

8.5 Comparison with Other Analyses 

Our results for e'/e cover the range of values between those obtained by the Dortmund 



17], [T8[ and Rome |19|] groups. Depending on the choice of the parameters listed in ( |3.12|) 



one can obtain the results of these two groups. 
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Here we would like to make a comparison with ref. |TS| because this is the only other 
paper in which full next-to-leading order corrections to the AS = 1 Hamiltonian have been 
computed. 

Since these authors presented their results only for the HV scheme, only the comparison 
for this scheme is possible. Here are our findings: 



i) The results for y 6 and y$ obtained in JU| agree with ours to better than 3%. Since the 
details on the two-loop anomalous dimensions are not yet available from the Rome 
group it is impossible for us to identify the origin of this difference^]. Fortunately 
this is immaterial for any phenomenological applications. We should however remark 
that the next-to-leading corrections in the HV scheme found by us are larger than 



those given in ref. [IjJ because there some of these corrections have been included in 
the leading order. 

ii) From the smooth //-dependence of y§ presented in jn| we conclude that the matching 
effects discussed in section 4.3 have not been taken into account there. In the case 
of the HV scheme as seen in fig. § this is fortunately only a small effect. 

Next we would like to list several differences between our own analysis and the one of 



ref. [ 1 9 1 which are relevant when taken individually but which cancel each other to some 



extent when taken as a whole. 

iii) Our analysis of section 6.5 shows that the contributions of Q\ 2 operators to e'/e 
evaluated at fi ~ 2 GeV lie at the lower end of the educated guess of ref. [ji~9] . We 
find Vt\ « -0.01, Vt c 2 « -0.04 to be compared with VL\ w 0.02, Vt c 2 « -0.18 quoted in 



[0|. This reduces the value of e'/e found by the Rome group. 



iv) Our analysis of B\ shows that the lattice calculations overestimate the contribu- 
tions of the operators Qg,io by almost a factor of two. When this is corrected for, the 
value of e' I e of ref. |H| is further reduced. 

v) The suppression of e'/e through the imposition of the AI =1/2 rule found by us is 



compatible with the one given in |L9| although somewhat stronger. 

3 ' We have just been informed by Guido Martinelli that the Rome group confirmed our NDR two- 
loop results of refs. JlJ, [is) . The final results on two-loop anomalous dimensions in the HV scheme also 
agree, although there are apparently some differences at the intermediate stages of the calculations. We 
guess then that the 3% difference in the Wilson coefficients originates in the numerical integration of the 
renormalization group equations done in [|19| which in our paper has been done analytically. 
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vi) As discussed in section 5 the values of m B (l GeV) used in |19| are higher than those 
used by us. With our values the results for e'/ e obtained in |19| would be increased. 



vii) The values of |Vub/Kb|, Bk and in particular of \V c b\ used in |19| are systematically 
higher than those used here which results in a smaller value of ImA t . With our choices 
of these parameters the authors of |B| would find a higher e'/e. In view of the 
increased 5-meson life time and the decrease of -Bk in the recent lattice calculations 



[pql we think that our estimate of ImA t is better than the one of ref. p9 |. 
As already stated above these effects cancel each other to some extent when taken 



as a whole and consequently our result given in tab. for e'/e in the HV scheme with 
Bq = B 8 = 1 are certainly consistent with the ones obtained in ref. [ITUJ. We have however 
varied the parameters Bq and B$ in a larger range and did in addition the same analysis in 
the NDR scheme. Consequently we reached somewhat different conclusions than ref. |19| 
regarding the two experimental results for e'/e given in ( |1 . 1| ) . 



9 Wilson Coefficients for AB = 1 Decays 

This section can be viewed as the generalization of the leading order analyses of Grinstein 
| 59fl and Kramer and Palmer beyond the leading order approximation. In addition, 
the coefficients of electroweak penguin operators will be calculated. 

We will focus on the AB = 1, AC = AU = part of the effective Hamiltonian which is 
of particular interest for the study of CP violation, i.e., decays to CP self-conjugate final 
states. 

At the tree level the effective Hamiltonian is simply given by 

n eff (AB = i) = % E E (9-i) 

q=u,c q i=d,s 

The cases q' = d or s can be treated separately and have the same coefficients C;(/x). 
Specifying to q' — d and using the unitarity of the CKM matrix we find for /i w 0(rrib) 
(6 = V* b V ld ) 

H eff (AB = 1) = ^ I £ c ( Zl (n) 0\ + z 2 (f,) 0\ 



(9.2) 



- it (5>(aO0<) J, 
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where Oj, i = 3, . . . , 10, are obtained from Q i} i = 3, . . . , 10 in ( |2.1[ ) by replacing s by 6, 

O" = (& a M/3) y _ A (M/3rf a ) y _ A ; O2 = (bu) v _ A (ud)y_ A (9.3) 



and the are obtained from (|9.3|) through the replacement u — > c. 

In spite of the presence of the CKM factors in ( |9.2| ), it is an easy matter to convince 
oneself that for /i = mj all ten coefficients building the vector C(/x) can be found using 

C(m b ) = Ur (m b ,M w ,a)C(M w ) , (9.4) 

with the evolution matrix of eq. Q2.17D given in an effective theory with 5 flavours and the 
initial conditions given in eqs. ([4.8|) - (|4.17 ). The numerical results for m t = 130 GeV are 



given in tab. |22|. They have a structure similar to the results in tab. |3| and therefore do 
not require further discussion. 

10 Ten Messages from This Paper 



Using the two-loop anomalous dimension matrices of refs. fT^- we have constructed 
the effective Hamiltonian for AS" = 1, AC = 1, and AB = 1 transitions with the Wilson 
coefficients Ci(fi) including leading and next-to- leading QCD corrections and leading order 
corrections in the electromagnetic coupling constant a. As the main application of the 
constructed AS = 1 Hamiltonian, we have calculated the ratio e'/e as a function of m t , 
Aj^g and of a set of non-perturbative parameters Bi. This new anatomy of e'/e improves 
the leading order anatomy of 1989 || in two respects: 

• The Wilson coefficients Cj(/i) include next-to- leading logarithmic corrections and 
consequently, the scale Ayg? extracted from deep-inelastic scattering and from e + e~- 
data at LEP can be meaningfully used for the first time in the evaluation of e'/ e. Since 
e' I e, being dominated by penguin contributions with large anomalous dimensions, is 
sensitive to the strength of the QCD coupling constant, the meaningful use of A^ 
increases considerably the precision of the theory provided the QCD scale can be 
extracted accurately from high energy data. 

• The hadronic matrix elements (Qi(yu)) used in the analysis are at all stages consis- 
tent with the data on CP-conserving K — » ttti decays. In particular, AI = 1/2 
enhancement and AI = 3/2 suppression are exactly built into our analysis. 

We have presented a detailed renormalization group analysis of both the Wilson coef- 
ficient functions and of the hadronic matrix elements. In particular, we have stressed that 
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the Wilson coefficients are renormalization scheme dependent and that this scheme depen- 
dence can only be cancelled by the one present in (Qi (//)). Since the existing methods 
are not capable to study this scheme and scale dependence accurately enough, we have 
presented a semi-phenomenological approach to hadronic matrix elements with the help of 
which we could determine several of these matrix elements from existing CP-conserving 
data. 

To summarize, the ten main messages of our paper are as follows: 

i) The coefficients 2_(AJ = 1/2) and z + (AI = 3/2) calculated in the NDR and HV 
schemes are contrary to the common wisdom suppressed and enhanced by next-to- 
leading order corrections, respectively. These effects must be compensated by those 
present in the hadronic matrix elements so that the experimentally known AI =1/2 
and AI = 3/2 amplitudes are reproduced. 

ii) The contribution of QCD penguins to the A J = 1/2 rule depends strongly on fi and 
the renormalization scheme considered, but is rather small with respect to the whole 
amplitude. For /i > m c it vanishes, but even for /i « 1 GeV it does not exceed 20% 
of A for m s (l GeV) > 125 MeV. In view of these features, relating the size of e'/e 
to the penguin contribution to A is in our opinion not a good idea. 

iii) The coefficient yQ(m c ) is suppressed by roughly 15% in the HV scheme but it is 
enhanced for fx > m c by 10% in the NDR scheme relative to the leading order value. 
Simulaneously, for m t < 150 GeV we find a strong enhancement of 2/s( m c) calculated 
in the HV scheme, to be compared with modest enhancement of y8{m c ) obtained in 
the NDR scheme. These enhancements of y$ become smaller for large m t values. 

iv) The choice fx = m c allows to extract efficiently several hadronic matrix elements from 
the data on CP-conserving K — > -kit decays. Indeed, at this scale the contributions 
of all penguin operators and furthermore Q\ 2 operators to the amplitudes A and 
A 2 vanish for the HV scheme or are very small for the NDR scheme. This allows 
for obtaining the matrix elements (Q 1,2(1^0)) 2 directly from A 2 as functions of A^jg. 
The matrix elements (<5i,2(^c))o are given as a function of B^ 2 \m c ) which is found 
to be in the range 5 < B^ 2 \m c ) < 7. The relations between operators allow to 
obtain subsequently all matrix elements of (V — A) <g> (V — A) penguin operators. 
This analysis shows that the extracted matrix elements differ considerably from the 
ones obtained by the vacuum insertion method and in the case of (Qi)2, i — 1,2, 9, 10 
also from the lattice results. We have also demonstrated that (Qi 2 (/•*)) 0,2 can be 
calculated and that for fx < 2 GeV the role of these operators in e'/e is small. 
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v) Our renormalization group analysis of (Qi{n)) shows an interesting result that the Bi 
parameters of all (V — A) <g> (V + A) operators are almost /i-independent. All P; pa- 
rameters of {V — A)®{V — A) operators show a substantial /i-dependence. The \i- and 
scheme dependences of (<5i,2 (/•*)) 2 and (Qg.io (/•*)) 2 are given to a very good approxima- 
tion by 1/z+(jj,) in the range 1 GeV < \i < 2 GeV. The /i-dependences of (Qs(ju))o, 
(Q6(A i ))o) (Q 7(^)2 and (<5s(a0)2 are given within 5% accuracy by l/m 2 (/V). We also 

(1/2) (3/2) 

point out that the contributions P 6 (/1) and P 8 (/1), which are /i-independent 
in the large-N limit show visible /^-dependences when exact renormalization group 
analysis with mixing between operators is taken into account. 

vi) As a preparation for future developments in the evaluation of hadronic matrix el- 
ements, we have presented e'/e in form of an expansion in the independent non- 
perturbative parameters Bi(m c ). The coefficients of this expansion are functions of 
renormalization scheme, Aj^g and m t and the values given in tabs. and |15| are such 



that the CP-conserving A and A 2 amplitudes are always consistent with the data. 

vii) Among several contributions to e'/e there are four which govern this ratio. Two of 
them can to a good accuracy be fixed by the CP-conserving data and consequently 
are only weakly dependent on A^ and the renormalization scheme. These are the 
negative m t -independent contribution of Q4 and a smaller positive contribution of 
Q9 + Q10 which moderately increases with m t . The former contribution is substan- 
tially enhanced when the AI = 1/2 rule is taken into account as pointed out already 
in ref. @, and recently reemphasized by the authors of ||19|| . On the other hand, Qg 
and Q10 contributions are similarly to A 2 suppressed substantially relatively to their 
vacuum insertion and lattice estimates. Consequently, if only the contributions of 
(V — A) ® (V — A) penguins Q4, Qg and Q±o are taken into account e'/e is negative 
and of O(10- 4 ). 

viii) The main uncertainties in e'/e result from ImA t and from the contributions of the 
{V — A) g) (V + A) penguin operators Q 6 and Q 8 which enter e'/e with opposite signs 
and are sensitive functions of A^jg. None of them can be determined from the existing 
CP-conserving data. Their matrix elements depend on m s and the parameters Pg 1 ^ 
and Pf /2) . For m s (m c ) = 150 MeV and P^ 1/2) (m c ) = Pf /2) (m c ) = 1, the Q 6 
contribution dominates for m t < 200 GeV, although this dominance is substantially 
suppressed in the HV scheme and moderately in the NDR scheme (see message iii) 
above). For m t < 150 GeV the contribution of Q$ is smaller than the one of C/4 
and consequently, for m t ~ (9(130 GeV) the ratio e'/e is mainly determined by the 
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QCD penguins Qq and Q4 with Q 6 being stronger and giving a positive e'/e. With 
increasing m t the role of Q 8 increases rapidly and for m t m (9(200 GeV) it is the 
competition of Qq and Q 8 which dominates e'/e. Q 8 being helped by the negative 
contribution of Q4 wins this competition for m t > 200 GeV and e' / e becomes negative 
similarly to the results of 0. 

Although the parameters B 6 and B 8 are essentially /^-independent and in accordance 
with the vacuum insertion method and the 1/N approach, there is an important 
question whether their actual values are in the neighborhood of unity. Recall that 
other parameters Bi show substantial deviation from Bi = 1. In this respect, we 



sympathize with the efforts made by the authors of ref. [17], y| [54J] who questioned 
the use of B 6 = B 8 = 1. Yet, in our opinion the chapter on the actual values of 
B 6 and B 8 is just being opened and it may well turn out that these parameters are 
indeed special with values close to unity as given in the leading order of the 1/N 
approach, in the vacuum insertion method, and in lattice calculations. A closely 
related issue is the actual value of m s . 

Our numerical analysis of e' / e summarized in tabs. |16|-pl| shows that for B e = B 8 = 1 
and high values of m t , at which strong cancellation takes place, the actual values for 
e'/e depend sensitively on whether the ratio y 8 /y§ is calculated in the NDR or HV 
scheme. This demonstrates very clearly the necessity of a solid evaluation of (Qe)o 
and (Q 8 )2 which would cancel this dependence. For higher values of B 6 and lower 
values of m t this scheme dependence is weaker. Generally, the values of e'/ e obtained 
in the HV scheme are below the ones obtained in the NDR scheme and roughly in 



agreement with the results of ref. fl9fl . We have also stressed that e'/e depends 
sensitively on A^. 

In the NDR scheme for B^ /2 \m c ) = B { 8 /2 \m c ) = 1, m s (m c ) = 150 MeV, A^ = 
300 MeV, and m t = 130 GeV, we find 

I (6.7 ± 2.6) x 10- 0<S< h 
1 \ (4.8 ±2.2) x 10~ 4 f < 5 < 7T, 

(1/2) (3/2) 

in perfect agreement with the E731 result. However, for Bq (m c ) = 2 and B 8 (m c ) 
1 we find 

f (20.0 ±6.5) xlO- 4 0<5<|, 
\ (14.4 ±5.6) x 10~ 4 !<<5<7r, 
in agreement with the findings of NA31. The LO results are slightly higher than 
these values. In the HV scheme somewhat lower values are found. 
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The fate of theoretical estimates of e'/e in the coming years depends crucially on 
whether it will be possible to reduce the uncertainties in the values of m t , Aj^g, m s , B^ 2 \ 
B^ 2 \ B K , | V^b | and | Kb/Kb I- If the standard CP violating interactions parametrized by 
the CKM matrix are the main origin of e' 7^ also the fate of experimental searchers for 
e'/e' 7^ depends crucially on these parameters. They should pray that the top quark will 
be found soon or that B 6 > 1.5, m s (l GeV) < 175 MeV and A m > 250 MeV. 
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Appendices 

A Quark Threshold Matching Matrices 



The matrix Sf s in ([4.36 ) has the following non- vanishing rows. 
(jl = m b : 

5r s (Q 4 ) = 5r s (Q 6 ) = -2 5f s (Q 8 ) = -2 5f s (Q 10 ) = -\ P (A.l) 

/I = m c : 

5r s (Q 4 ) = 5r s (Q 6 ) = 8r 9 (Qs) = Sr s (Qi ) = ~P (A.2) 

y 

where 

P = (0,0,-^,1, --,1,0,0,0,0). (A.3) 

The matrix 5r e has the following non- vanishing rows. 
\i = mb: 

5r c (Q 3 ) = 3 5r e (Qd = Sr e (Qs) = 3 5r c (Q e ) = ^ P (A.4) 



5f e (Q 7 ) = 3 5f c (Q 8 ) = 5f c (Q 9 ) = 3 5f c (Q w ) = ~P (A.5) 



27 
1( 
27 



// = m c : 

Sf c (Q 3 ) = Sr c (Q 5 ) = Sf c (Q 7 ) = 5f c (Q 9 ) = ~P (A.6) 



27 

5r e (Q 4 ) = *r e (Q 6 ) = 5f c (Q 8 ) = 5f c (Q 10 ) = ~P (A.7) 

ol 



where 

P= (0,0,0,0,0,0,1,0,1,0). (A.i 
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B Numerical Results for Wilson Coefficients, e 1 je and 
Bi- Expansion 

Table 1: AS = 1 Wilson coefficients at /j, — 1 GeV for m t = 130 GeV. In our approach 
y x — y 2 = holds. 





A M s = 0.2GeV 


Aj^t = 0.3 GeV 


= 0.4 GeV 


Scheme 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


z\ 


-0.587 


-0.397 


-0.477 


-0.715 


-0.486 


-0.606 


-0.854 


-0.588 


-0.774 


Z2 


1.319 


1.204 


1.256 


1.409 


1.262 


1.345 


1.511 


1.333 


1.472 


Z3 


0.004 


0.008 


0.004 


0.005 


0.013 


0.008 


0.007 


0.021 


0.015 


Zi 


-0.010 


-0.023 


-0.011 


-0.014 


-0.035 


-0.018 


-0.018 


-0.054 


-0.030 


^5 


0.003 


0.006 


0.003 


0.004 


0.007 


0.004 


0.006 


0.009 


0.006 


Z(y 


-0.011 


-0.023 


-0.010 


-0.015 


-0.035 


-0.016 


-0.020 


-0.056 


-0.026 


z-j/a 


0.005 


0.000 


-0.005 


0.008 


0.008 


-0.003 


0.012 


0.017 


-0.002 


z 8 /a 


0.001 


0.009 


0.007 


0.001 


0.016 


0.011 


0.002 


0.028 


0.018 


z 9 /a 


0.005 


0.005 


-0.001 


0.009 


0.016 


0.004 


0.013 


0.030 


0.010 


zio/a 


-0.001 


-0.006 


-0.007 


-0.001 


-0.009 


-0.011 


-0.001 


-0.015 


-0.017 


Z- 


1.907 


1.600 


1.733 


2.123 


1.748 


1.951 


2.365 


1.921 


2.245 


z+ 


0.732 


0.807 


0.779 


0.694 


0.776 


0.740 


0.657 


0.746 


0.698 


z-/z+ 


2.604 


1.982 


2.226 


3.060 


2.252 


2.638 


3.597 


2.576 


3.217 




ys 


0.027 


0.022 


0.025 


0.034 


0.029 


0.033 


0.042 


0.036 


0.043 


y4 


-0.048 


-0.044 


-0.046 


-0.056 


-0.052 


-0.055 


-0.064 


-0.060 


-0.065 


y<o 


0.011 


0.005 


0.013 


0.012 


0.001 


0.015 


0.013 


-0.008 


0.018 


y% 


-0.078 


-0.071 


-0.065 


-0.102 


-0.099 


-0.087 


-0.130 


-0.142 


-0.118 


yi/a 


-0.025 


-0.083 


-0.082 


-0.017 


-0.080 


-0.079 


-0.009 


-0.078 


-0.076 


y$/a 


0.053 


0.069 


0.081 


0.074 


0.095 


0.110 


0.101 


0.133 


0.154 


y 9 /a 


-1.160 


-1.161 


-1.161 


-1.222 


-1.225 


-1.226 


-1.297 


-1.314 


-1.318 


yio/a 


0.488 


0.400 


0.407 


0.592 


0.502 


0.512 


0.706 


0.634 


0.647 
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Table 2: AS = 1 Wilson coefficients at /i = 1.4 GeV for m t = 130 GeV. In our approach 
y x — y 2 = holds. 





A^ = 0.2GeV 


A^ = 0.3 GeV 


= 0.4 GeV 


Scheme 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 




-0.486 


-0.321 


-0.383 


-0.573 


-0.379 


-0.459 


-0.660 


-0.438 


-0.543 


^2 


1.252 


1.158 


1.195 


1.310 


1.193 


1.244 


1.369 


1.231 


1.301 


Z- 


1.738 


1.479 


1.578 


1.883 


1.571 


1.704 


2.030 


1.669 


1.844 


~r 


0.766 


0.836 


0.812 


0.736 


0.814 


0.785 


0.709 


0.793 


0.758 




2.269 


1.769 


1.943 


2.558 


1.930 


2.170 


2.862 


2.105 


2.432 




2/3 


0.024 


0.020 


0.022 


0.029 


0.025 


0.028 


0.035 


0.030 


0.035 


2/4 


-0.047 


-0.044 


-0.045 


-0.055 


-0.052 


-0.054 


-0.063 


-0.060 


-0.063 


2/5 


0.013 


0.007 


0.013 


0.014 


0.006 


0.015 


0.015 


0.004 


0.017 


2/6 


-0.068 


-0.061 


-0.056 


-0.085 


-0.077 


-0.070 


-0.102 


-0.097 


-0.087 


2/7/a 


-0.023 


-0.082 


-0.081 


-0.014 


-0.080 


-0.079 


-0.006 


-0.078 


-0.076 


2/s/a 


0.041 


0.059 


0.067 


0.054 


0.076 


0.085 


0.068 


0.097 


0.106 


2/g/a 


-1.103 


-1.110 


-1.110 


-1.139 


-1.144 


-1.143 


-1.179 


-1.182 


-1.182 


2/io/a 


0.403 


0.320 


0.325 


0.474 


0.379 


0.385 


0.545 


0.443 


0.450 



Table 3: AS = 1 Wilson coefficients at /i = 2 GeV for m t = 130 GeV. In our approach 
y l — y 2 = holds. 





A m = 0.2 GeV 


A^ = 0.3 GeV 


A^ = 0.4 GeV 


Scheme 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


Zl 


-0.403 


-0.261 


-0.311 


-0.465 


-0.300 


-0.362 


-0.523 


-0.338 


-0.413 


Z2 


1.200 


1.123 


1.152 


1.238 


1.145 


1.182 


1.276 


1.168 


1.214 


Z- 


1.603 


1.384 


1.463 


1.703 


1.445 


1.545 


1.798 


1.505 


1.627 


z+ 


0.797 


0.862 


0.840 


0.773 


0.845 


0.820 


0.752 


0.830 


0.801 


Z-/Z+ 


2.010 


1.606 


1.741 


2.202 


1.710 


1.884 


2.390 


1.814 


2.030 




2/3 


0.019 


0.019 


0.018 


0.023 


0.022 


0.021 


0.026 


0.026 


0.025 


2/4 


-0.039 


-0.045 


-0.038 


-0.045 


-0.052 


-0.044 


-0.050 


-0.060 


-0.050 


2/5 


0.011 


0.010 


0.011 


0.012 


0.010 


0.013 


0.013 


0.011 


0.014 


L 2/6 


-0.053 


-0.056 


-0.045 


-0.064 


-0.067 


-0.054 


-0.074 


-0.080 


-0.063 


2/7/a 


-0.023 


-0.077 


-0.073 


-0.016 


-0.073 


-0.072 


-0.010 


-0.069 


-0.071 


2/s/a 


0.032 


0.047 


0.054 


0.040 


0.057 


0.065 


0.048 


0.068 


0.076 


y 9 /a 


-1.060 


-1.073 


-1.070 


-1.084 


-1.092 


-1.091 


-1.108 


-1.111 


-1.114 


yw/a 


0.333 


0.264 


0.268 


0.384 


0.304 


0.309 


0.432 


0.344 


0.349 
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Table 4: AS = 1 Wilson coefficient y-jja at ji — 1 GeV for various m t values. 





A m = 0.2 GeV 


A m = 0.3 GeV 


= 0.4 GeV 


m t [ GeV] 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


110 


-0.045 


-0.102 


-0.102 


-0.036 


-0.100 


-0.099 


-0.028 


-0.097 


-0.094 


130 


-0.025 


-0.083 


-0.082 


-0.017 


-0.080 


-0.079 


-0.009 


-0.078 


-0.076 


150 


-0.002 


-0.058 


-0.058 


0.006 


-0.057 


-0.056 


0.014 


-0.056 


-0.053 


170 


0.025 


-0.031 


-0.031 


0.033 


-0.030 


-0.029 


0.039 


-0.030 


-0.027 


190 


0.056 


-0.000 


-0.000 


0.062 


-0.000 


0.001 


0.068 


-0.001 


0.002 



Table 5: AS = 1 Wilson coefficient y%ja at /i — 1 GeV for various m t values. 





A m = 0.2 GeV 


A^ = 0.3 GeV 


A^=0.4 GeV 


m t [GeV] 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


110 


0.031 


0.050 


0.062 


0.045 


0.070 


0.085 


0.062 


0.098 


0.119 


130 


0.053 


0.069 


0.081 


0.074 


0.095 


0.110 


0.101 


0.133 


0.154 


150 


0.079 


0.090 


0.102 


0.109 


0.124 


0.140 


0.145 


0.174 


0.195 


170 


0.109 


0.115 


0.126 


0.148 


0.158 


0.173 


0.196 


0.221 


0.242 


190 


0.142 


0.142 


0.153 


0.191 


0.195 


0.210 


0.252 


0.273 


0.293 



Table 6: AS = 1 Wilson coefficient sum (y g + yio)/a at ji — 1 GeV for various m t values. 





A^ = 0.2 GeV 


A^ = 0.3 GeV 


A M g = 0.4 GeV 


m t [GeV] 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


110 


-0.567 


-0.652 


-0.646 


-0.531 


-0.620 


-0.613 


-0.497 


-0.585 


-0.576 


130 


-0.672 


-0.760 


-0.754 


-0.630 


-0.722 


-0.715 


-0.591 


-0.680 


-0.671 


150 


-0.779 


-0.872 


-0.866 


-0.732 


-0.828 


-0.820 


-0.688 


-0.779 


-0.770 


170 


-0.891 


-0.987 


-0.981 


-0.838 


-0.937 


-0.929 


-0.788 


-0.881 


-0.872 


190 


-1.008 


-1.108 


-1.102 


-0.948 


-1.051 


-1.043 


-0.893 


-0.988 


-0.978 
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Table 7: AS = 1 Wilson coefficient 1/7/ a at /i — 1.4 GeV for various m t values. 





= 0.2 GeV 


A^ = 0.3 GeV 


A^ = 0.4 GeV 


m t [GeV] 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


110 


-0.043 


-0.102 


-0.102 


-0.034 


-0.100 


-0.099 


-0.025 


-0.098 


-0.096 


130 


-0.023 


-0.082 


-0.081 


-0.014 


-0.080 


-0.079 


-0.006 


-0.078 


-0.076 


150 


0.001 


-0.057 


-0.057 


0.010 


-0.056 


-0.054 


0.018 


-0.054 


-0.052 


170 


0.029 


-0.029 


-0.029 


0.037 


-0.028 


-0.027 


0.044 


-0.027 


-0.025 


190 


0.060 


0.002 


0.003 


0.067 


0.003 


0.004 


0.074 


0.003 


0.005 



Table 8: AS = 1 Wilson coefficient y^/a. at ji — 1.4 GeV for various m t values. 





A M 3 = 0.2GeV 


Aj^t = 0.3 GeV 


A M ^=0.4 GeV 


m t [GeV] 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


110 


0.024 


0.045 


0.053 


0.032 


0.059 


0.068 


0.042 


0.075 


0.085 


130 


0.041 


0.059 


0.067 


0.054 


0.076 


0.085 


0.068 


0.097 


0.106 


150 


0.061 


0.075 


0.084 


0.079 


0.097 


0.106 


0.099 


0.122 


0.132 


170 


0.084 


0.094 


0.102 


0.107 


0.120 


0.129 


0.133 


0.151 


0.161 


190 


0.109 


0.115 


0.123 


0.139 


0.146 


0.155 


0.172 


0.183 


0.193 



Table 9: AS = 1 Wilson coefficient sum (yg + yio)/ct at /i — 1.4 GeV for various m t values. 





= 0.2 GeV 


A^ = 0.3 GeV 


A^ = 0.4 GeV 


m t [GeV] 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


110 


-0.591 


-0.677 


-0.672 


-0.560 


-0.655 


-0.649 


-0.532 


-0.634 


-0.627 


130 


-0.700 


-0.791 


-0.785 


-0.665 


-0.764 


-0.758 


-0.633 


-0.739 


-0.732 


150 


-0.812 


-0.907 


-0.902 


-0.773 


-0.877 


-0.871 


-0.738 


-0.847 


-0.840 


170 


-0.929 


-1.028 


-1.023 


-0.886 


-0.993 


-0.987 


-0.846 


-0.959 


-0.952 


190 


-1.052 


-1.155 


-1.149 


-1.003 


-1.115 


-1.109 


-0.959 


-1.076 


-1.069 
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Table 10: AS = 1 Wilson coefficient y-jja at /i — 2 GeV for various m t values. 





= 0.2 GeV 


A^ = 0.3 GeV 


A^ = 0.4 GeV 


m t [GeV] 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


110 


-0.043 


-0.098 


-0.094 


-0.036 


-0.093 


-0.093 


-0.030 


-0.089 


-0.091 


130 


-0.023 


-0.077 


-0.073 


-0.016 


-0.073 


-0.072 


-0.010 


-0.069 


-0.071 


150 


0.002 


-0.051 


-0.048 


0.008 


-0.048 


-0.047 


0.014 


-0.044 


-0.046 


170 


0.030 


-0.023 


-0.020 


0.036 


-0.020 


-0.019 


0.041 


-0.017 


-0.018 


190 


0.062 


0.009 


0.013 


0.067 


0.012 


0.013 


0.072 


0.015 


0.013 



Table 11: AS = 1 Wilson coefficient ys/at at fi = 2 GeV for various m t values. 





A M 3 = 0.2GeV 


Aj^t = 0.3 GeV 


A M ^=0.4 GeV 


m t [GeV] 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


110 


0.019 


0.036 


0.043 


0.024 


0.044 


0.052 


0.029 


0.052 


0.061 


130 


0.032 


0.047 


0.054 


0.040 


0.057 


0.065 


0.048 


0.068 


0.076 


150 


0.048 


0.059 


0.067 


0.059 


0.072 


0.080 


0.070 


0.086 


0.095 


170 


0.066 


0.074 


0.081 


0.081 


0.090 


0.098 


0.096 


0.106 


0.115 


190 


0.086 


0.090 


0.098 


0.104 


0.109 


0.117 


0.124 


0.129 


0.138 



Table 12: AS = 1 Wilson coefficient sum (y g +yio)/at at /i — 2 GeV for various m t values. 





= 0.2 GeV 


A^ = 0.3 GeV 


A^ = 0.4 GeV 


m t [GeV] 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


110 


-0.613 


-0.693 


-0.685 


-0.589 


-0.674 


-0.669 


-0.569 


-0.656 


-0.653 


130 


-0.727 


-0.810 


-0.802 


-0.700 


-0.788 


-0.782 


-0.676 


-0.767 


-0.764 


150 


-0.844 


-0.931 


-0.923 


-0.814 


-0.905 


-0.900 


-0.787 


-0.882 


-0.879 


170 


-0.966 


-1.056 


-1.048 


-0.932 


-1.028 


-1.022 


-0.902 


-1.001 


-0.998 


190 


-1.093 


-1.187 


-1.179 


-1.055 


-1.155 


-1.150 


-1.022 


-1.125 


-1.123 
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Table 13: AS = 1 Wilson coefficient ratio dml2±mlll a t // = 1.4 GeV for various m t 
values. 





0.2 GeV 


A M ^ = 0.3 GeV 


A^=0.4 GeV 


m t [GeV] 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


110 


0.144 


0.183 


0.349 


0.252 


0.331 


0.499 


0.330 


0.443 


0.615 


130 


0.492 


0.525 


0.714 


0.582 


0.644 


0.838 


0.646 


0.730 


0.935 


150 


0.898 


0.923 


1.138 


0.964 


1.007 


1.233 


1.013 


1.066 


1.308 


170 


1.352 


1.370 


1.619 


1.395 


1.418 


1.680 


1.426 


1.445 


1.729 


190 


1.855 


1.867 


2.146 


1.870 


1.871 


2.173 


1.883 


1.864 


2.195 



Table 14: Coefficients in the expansion of pW 2 \ 







LO 


NDR 


HV 




mt[ GeV] 


a (1/2) 
"0 


a (1/2) 


a (1/2) 
"6 


a (1/2) 




"6 


a (1/2) 
"0 


a 2 


"6 




110 


-6.49 


0.48 


9.84 


-7.66 


0.57 


8.97 


-7.59 


0.60 


8.03 




130 


-6.18 


0.45 


9.97 


-7.29 


0.53 


9.07 


-7.22 


0.56 


8.13 


200 


150 


-5.84 


0.41 


10.07 


-6.89 


0.49 


9.15 


-6.82 


0.52 


8.20 




170 


-5.46 


0.38 


10.15 


-6.45 


0.45 


9.21 


-6.38 


0.48 


8.26 




190 


-5.05 


0.34 


10.21 


-5.97 


0.41 


9.26 


-5.90 


0.43 


8.32 




110 


-6.62 


0.47 


12.29 


-7.94 


0.57 


11.56 


-7.84 


0.61 


10.07 




130 


-6.33 


0.44 


12.44 


-7.61 


0.54 


11.68 


-7.51 


0.57 


10.19 


300 


150 


-6.00 


0.42 


12.56 


-7.24 


0.51 


11.77 


-7.14 


0.54 


10.28 




170 


-5.65 


0.39 


12.66 


-6.84 


0.47 


11.85 


-6.74 


0.50 


10.36 




190 


-5.26 


0.35 


12.74 


-6.41 


0.43 


11.91 


-6.31 


0.46 


10.42 




110 


-6.68 


0.46 


14.94 


-8.13 


0.57 


14.78 


-8.00 


0.60 


12.49 




130 


-6.40 


0.44 


15.13 


-7.82 


0.54 


14.92 


-7.69 


0.57 


12.64 


400 


150 


-6.08 


0.41 


15.27 


-7.48 


0.51 


15.03 


-7.35 


0.54 


12.75 




170 


-5.73 


0.39 


15.39 


-7.10 


0.48 


15.12 


-6.97 


0.51 


12.84 




190 


-5.35 


0.36 


15.48 


-6.70 


0.45 


15.19 


-6.57 


0.48 


12.91 
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Table 15: Coefficients in the expansion of P^ 2 \ 







LO 


NDR 


HV 


^MS 


m t [GeV] 


fl (3/2) 

"o 


a (3/2) 
"8 




a (3/2) 


a (3/2) 


a (3/2) 




110 


-1.01 


0.66 


-1.04 


0.87 


-1.08 


1.35 




130 


-1.19 


1.97 


-1.22 


2.01 


-1.26 


2.49 


200 


150 


-1.38 


3.52 


-1.41 


3.36 


-1.45 


3.84 




170 


-1.58 


5.28 


-1.60 


4.89 


-1.64 


5.37 




190 


-1.79 


7.24 


-1.80 


6.59 


-1.84 


7.07 




110 


-0.99 


1.29 


-1.03 


1.68 


-1.08 


2.22 




130 


-1.18 


2.84 


-1.21 


3.02 


-1.26 


3.55 


300 


150 


-1.37 


4.67 


-1.40 


4.59 


-1.45 


5.13 




170 


-1.57 


6.76 


-1.59 


6.38 


-1.64 


6.92 




190 


-1.78 


9.07 


-1.79 


8.37 


-1.84 


8.91 




110 


-0.98 


1.96 


-1.02 


2.66 


-1.08 


3.24 




130 


-1.16 


3.78 


-1.20 


4.23 


-1.26 


4.81 


400 


150 


-1.36 


5.92 


-1.39 


6.07 


-1.44 


6.66 




170 


-1.56 


8.36 


-1.58 


8.17 


-1.64 


8.75 




190 


-1.76 


11.06 


-1.78 


10.50 


-1.84 


11.09 
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Table 16: e'/e in units of 1(T 4 for B ( q /2) = 1.0 and B£ /2) = 1.0. 







LO 


NDR 


HV 




m t 


1. Quad. 


2. Quad. 


1. Quad. 


2. Quad. 


1. Quad. 


2. Quad. 




130 


6.4 ± 2.4 


4.6 ± 2.0 


5.2 ± 2.2 


3.7 ± 1.8 


3.7 ± 1.7 


2.6 ± 1.4 


200 


150 


4.6 ± 1.9 


3.2 ± 1.6 


3.7 ± 1.7 


2.6 ± 1.4 


2.4 ± 1.3 


1.6 ± 1.0 




170 


3.0 ± 1.3 


2.0 ± 1.1 


2.3 ± 1.2 


1.6 ± 1.0 


1.1 ± 0.8 


0.8 ± 0.6 




130 


8.0 ± 2.9 


5.8 ± 2.5 


6.7 ± 2.6 


4.8 ± 2.2 


4.6 ± 2.0 


3.3 ± 1.6 


300 


150 


5.8 ± 2.2 


4.1 ± 1.9 


4.9 ± 2.0 


3.4 ± 1.7 


2.9 ± 1.5 


2.0 ± 1.2 




170 


3.8 ± 1.6 


2.6 ± 1.4 


3.1 ± 1.5 


2.1 ± 1.2 


1.4 ± 0.9 


1.0 ± 0.7 




130 


9.9 ± 3.4 


7.1 ± 2.9 


8.8 ± 3.2 


6.3 ± 2.7 


5.7 ± 2.3 


4.1 ± 1.9 


400 


150 


7.2 ± 2.7 


5.0 ± 2.3 


6.5 ± 2.5 


4.5 ± 2.2 


3.7 ± 1.7 


2.6 ± 1.4 




170 


4.7 ± 1.9 


3.2 ± 1.6 


4.3 ± 1.9 


3.0 ± 1.6 


1.8 ± 1.1 


1.3 ± 0.9 



Table 17: e'/e in units of 10" 4 for Eg m = 1.5 and Bf ,2) = 1.5. 







LO 


NDR 


HV 




mt 


1. Quad. 


2. Quad. 


1. Quad. 


2. Quad. 


1. Quad. 


2. Quad. 




130 


10.9 ± 3.8 


7.9 ± 3.2 


9.2 ± 3.3 


6.6 ± 2.8 


6.9 ± 2.7 


5.0 ± 2.2 


200 


150 


8.0 ± 2.9 


5.6 ± 2.5 


6.6 ± 2.6 


4.6 ± 2.2 


4.6 ± 1.9 


3.2 ± 1.6 




170 


5.2 ± 2.1 


3.6 ± 1.8 


4.3 ± 1.8 


2.9 ± 1.6 


2.5 ± 1.2 


1.7 ± 1.0 




130 


13.5 ± 4.5 


9.7 ± 3.9 


11.6 ± 4.1 


8.3 ± 3.4 


8.3 ± 3.1 


6.0 ± 2.6 


300 


150 


9.9 ± 3.5 


6.9 ± 3.0 


8.5 ± 3.2 


5.9 ± 2.7 


5.6 ± 2.3 


3.9 ± 1.9 




170 


6.5 ± 2.5 


4.5 ± 2.2 


5.6 ± 2.3 


3.8 ± 2.0 


3.0 ± 1.4 


2.0 ± 1.2 




130 


16.3 ± 5.4 


11.8 ± 4.6 


14.8 ± 5.0 


10.6 ± 4.3 


10.2 ± 3.7 


7.3 ± 3.1 


400 


150 


12.0 ± 4.1 


8.4 ± 3.6 


11.0 ± 3.9 


7.7 ± 3.4 


6.9 ± 2.7 


4.8 ± 2.3 




170 


8.0 ± 3.0 


5.5 ± 2.6 


7.5 ± 2.9 


5.1 ± 2.5 


3.7 ± 1.7 


2.5 ± 1.4 
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Table 18: e'/e in units of 10~ 4 for b£ ,2) = 0.75 and Bf /2) = 0.75. 







LO 


NDR 


HV 




m t 


1. Quad. 


2. Quad. 


1. Quad. 


2. Quad. 


1. Quad. 


2. Quad. 




130 


4.1 ± 1.8 


3.0 ± 1.4 


3.2 ± 1.6 


2.3 ± 1.3 


2.1 ± 1.3 


1.5 ± 1.0 


200 


150 


2.9 ± 1.4 


2.0 ± 1.1 


2.2 ± 1.2 


1.6 ± 1.0 


1.2 ± 0.9 


0.9 ± 0.7 




170 


1.8 ± 1.0 


1.2 ± 0.8 


1.3 ± 0.9 


0.9 ± 0.7 


0.4 ± 0.6 


0.3 ± 0.4 




130 


5.3 ± 2.1 


3.8 ± 1.7 


4.3 ± 1.9 


3.1 ± 1.6 


2.7 ± 1.4 


1.9 ± 1.2 


300 


150 


3.8 ± 1.6 


2.6 ± 1.4 


3.1 ± 1.5 


2.1 ± 1.2 


1.6 ± 1.0 


1.1 ± 0.8 




170 


2.4 ± 1.2 


1.6 ± 1.0 


1.9 ± 1.1 


1.3 ± 0.9 


0.6 ± 0.7 


0.4 ± 0.5 




130 


6.6 ± 2.5 


4.8 ± 2.1 


5.8 ± 2.4 


4.2 ± 2.0 


3.5 ± 1.7 


2.5 ± 1.4 


400 


150 


4.8 ± 1.9 


3.4 ± 1.6 


4.2 ± 1.9 


2.9 ± 1.5 


2.2 ± 1.2 


1.5 ± 1.0 




170 


3.1 ± 1.4 


2.1 ± 1.2 


2.8 ± 1.4 


1.9 ± 1.1 


0.9 ± 0.8 


0.6 ± 0.6 



Table 19: e'/e in units of 10" 4 for Bg /2) = 1.25 and Bf' i) = 0.75. 







LO 


NDR 


HV 




mt 


1. Quad. 


2. Quad. 


1. Quad. 


2. Quad. 


1. Quad. 


2. Quad. 




130 


9.8 ± 3.4 


7.0 ± 2.9 


8.3 ± 3.1 


6.0 ± 2.6 


6.7 ± 2.6 


4.8 ± 2.2 


200 


150 


8.1 ± 2.9 


5.6 ± 2.5 


6.9 ± 2.7 


4.8 ± 2.3 


5.4 ± 2.2 


3.8 ± 1.9 




170 


6.6 ± 2.5 


4.5 ± 2.2 


5.6 ± 2.3 


3.8 ± 1.9 


4.3 ± 1.8 


2.9 ± 1.5 




130 


12.4 ± 4.2 


8.9 ± 3.6 


10.9 ± 3.9 


7.8 ± 3.3 


8.5 ± 3.2 


6.1 ± 2.7 


300 


150 


10.3 ± 3.6 


7.1 ± 3.1 


9.1 ± 3.3 


6.3 ± 2.9 


6.9 ± 2.7 


4.8 ± 2.3 




170 


8.3 ± 3.1 


5.7 ± 2.7 


7.4 ± 2.9 


5.0 ± 2.5 


5.4 ± 2.2 


3.7 ± 1.9 




130 


15.3 ± 5.0 


11.0 ± 4.3 


14.2 ± 4.9 


10.2 ± 4.1 


10.7 ± 3.8 


7.7 ± 3.2 


400 


150 


12.7 ± 4.3 


8.8 ± 3.8 


11.9 ± 4.2 


8.3 ± 3.7 


8.7 ± 3.2 


6.1 ± 2.8 




170 


10.3 ± 3.7 


7.0 ± 3.2 


9.8 ± 3.6 


6.6 ± 3.2 


6.9 ± 2.7 


4.7 ± 2.3 
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Table 20: e'/e in units of 10" 4 for B ( q /2) = 1.5 and Bf ,2) = 1.0. 







LO 


NDR 


HV 




m t 


1. Quad. 


2. Quad. 


1. Quad. 


2. Quad. 


1. Quad. 


2. Quad. 




130 


12.1 ± 4.1 


8.7 ± 3.5 


10.3 ± 3.7 


7.4 ± 3.1 


8.3 ± 3.1 


6.0 ± 2.6 


200 


150 


9.8 ± 3.5 


6.8 ± 3.0 


8.4 ± 3.1 


5.8 ± 2.7 


6.6 ± 2.6 


4.6 ± 2.2 




170 


7.7 ± 2.9 


5.2 ± 2.5 


6.6 ± 2.6 


4.5 ± 2.2 


5.0 ± 2.1 


3.4 ± 1.7 




130 


15.1 ± 5.0 


10.9 ± 4.3 


13.3 ± 4.6 


9.6 ± 3.9 


10.4 ± 3.7 


7.5 ± 3.1 


300 


150 


12.3 ± 4.2 


8.6 ± 3.7 


10.9 ± 3.9 


7.6 ± 3.4 


8.2 ± 3.1 


5.7 ± 2.6 




170 


9.7 ± 3.5 


6.6 ± 3.1 


8.6 ± 3.3 


5.9 ± 2.8 


6.2 ± 2.5 


4.2 ± 2.1 




130 


18.5 ± 6.0 


13.3 ± 5.1 


17.2 ± 5.7 


12.4 ± 4.9 


12.9 ± 4.5 


9.3 ± 3.8 


400 


150 


15.1 ± 5.1 


10.5 ± 4.5 


14.2 ± 4.9 


9.9 ± 4.3 


10.3 ± 3.7 


7.2 ± 3.2 




170 


11.9 ± 4.3 


8.1 ± 3.7 


11.3 ± 4.2 


7.7 ± 3.6 


7.8 ± 3.0 


5.3 ± 2.6 



Table 21: e'/e in units of 10" 4 for B\ 1 ' = 2.0 and Bf ' = 1.0. 







LO 


NDR 


HV 




m t 


1. Quad. 


2. Quad. 


1. Quad. 


2. Quad. 


1. Quad. 


2. Quad. 




130 


17.8 ± 5.8 


12.8 ± 5.0 


15.5 ± 5.2 


11.1 ± 4.4 


12.9 ± 4.5 


9.3 ± 3.8 


200 


150 


15.0 ± 5.0 


10.4 ± 4.4 


13.1 ± 4.5 


9.1 ± 4.0 


10.8 ± 3.8 


7.5 ± 3.3 




170 


12.4 ± 4.4 


8.5 ± 3.9 


10.9 ± 4.0 


7.4 ± 3.5 


8.8 ± 3.3 


6.0 ± 2.9 




130 


22.2 ± 7.1 


16.0 ± 6.1 


20.0 ± 6.5 


14.4 ± 5.6 


16.2 ± 5.4 


11.6 ± 4.6 


300 


150 


18.8 ± 6.2 


13.1 ± 5.5 


16.9 ± 5.7 


11.8 ± 5.0 


13.5 ± 4.7 


9.4 ± 4.1 




170 


15.6 ± 5.4 


10.6 ± 4.8 


14.1 ± 5.1 


9.6 ± 4.4 


11.1 ± 4.1 


7.5 ± 3.5 




130 


27.1 ± 8.5 


19.5 ± 7.4 


25.7 ± 8.2 


18.5 ± 7.1 


20.1 ± 6.6 


14.5 ± 5.6 


400 


150 


22.9 ± 7.5 


16.0 ± 6.6 


21.8 ± 7.2 


15.2 ± 6.4 


16.8 ± 5.7 


11.7 ± 5.0 




170 


19.1 ± 6.6 


13.0 ± 5.8 


18.3 ± 6.4 


12.5 ± 5.6 


13.8 ± 5.0 


9.4 ± 4.3 



72 



Table 22: AB = 1 Wilson coefficients at // = for m t = 130 GeV. 





A^ = 0.2GeV 


A^r = 0.3 GeV 


A m = 0.4 GeV 


Scheme 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 




-0.255 


-0.152 


-0.187 


-0.286 


-0.170 


-0.210 


-0.312 


-0.185 


-0.230 




1.115 


1.066 


1.083 


1.131 


1.075 


1.095 


1.146 


1.082 


1.106 


yz 


0.011 


0.011 


0.010 


0.012 


0.013 


0.011 


0.014 


0.014 


0.013 




-0.025 


-0.029 


-0.024 
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C Collection of Numerical Input Parameters 



Quark Masses 

m t = 100 - 200 GeV m centrai = 130 GeV 

m b = 4.8 GeV m c = 1.4 GeV 

m d (m c ) = 8 MeV m s (m c ) = 150 MeV 

Scalar Meson Masses and Decay Constants 

m K = 498 MeV F K = 161 MeV 

m w = 135 MeV = 132 MeV 



QCD and Electroweak Parameters 

A m = 200 - 400 MeV A^| tral = 300 MeV 

a = 1/128 M w = 80.0 GeV 

CKM Elements 

|Ks| = 0.221 |Kd| = 0.9753 

\V ch \ = 0.043 ±0.004 \V uh /V ch \ = 0.10 ±0.03 (range I) 

\V ch \ = 0.043 ±0.002 \V uh /V ch \ = 0.10 ±0.01 (range II) 

K — > 7T7T Decays and K° — K° Mixing 

ReA = 3.33 x 10" 7 GeV ReA 2 = 1.50 x lO" 8 GeV 

u = 1/22.2 tt vv , = 0.25 

e = (2.258 ± 0.018) x 10~ 3 AM K = 3.5 x lO" 15 GeV 

B K = 0.65 ±0.15 (range I) B K = 0.65 ± 0.05 (range II) 

?7i = 0.85 r] 2 = 0.58 

rj 3 = 0.36 

Hadronic Matrix Elements 

Bif (m c ) = 5.8 ±1.1 <ndrK) = 6.7 ± 0.9 <^K) = 6.3 ±1.0 

B (im = B (i/2) = B (i/2) = B (i m = s (i/ 2 ) = B w) = B w) = x (central values) 
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Figures 
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Figure 4: 
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Figure 11: 
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Figure 14: 
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Figure 16: 
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Figure 18: 
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Figure 19: 
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Figure Captions 

Figure [l|: Current-current and penguin 1-loop diagrams in the full theory. The wavy 
lines labeled with "W" denote W ± bosons, the unlabeled wavy lines denote either a 
gluon or a photon. 

Figure [2|: Current-current and penguin 1-loop diagrams in the effective theory. The 
unlabeled wavy lines denote either a gluon or a photon. 

Figure [3|: z±(fi) for various schemes and Aj^g = 300 MeV. 

Figure [4]: yeifi) and ys(^)/a for various schemes and Aj^g = 300 MeV. 

Figure [5|: y 7 (m c )/a and y%{m c )/a as functions of m t for various schemes and Aj^g = 
300 MeV. 

Figure [6|: yg(m c )/a and yio(m c )/a as functions of m t for various schemes and A^jg = 
300 MeV. 

Figure [7]: B^l 2 \n) entering (Qi,2(/-0)2 extracted from Rev4 2 for different schemes, A^jg = 
300 MeV and different cases discussed in the text. 

Figure [8|: £>( 3//2 )(/x) entering (<5g,io(A t ))2 extracted from ReA 2 for different schemes, A^jg = 
300 MeV and different cases discussed in the text. 

Figure |^: B± (m c ) as a function of B^ (m c ) for A^jg = 300 MeV and various schemes. 

Figure [T0| : B^ 2 ' as a function of B^ (m c ) for various values of B^ 2 \m c ). 



Figure [TT| : B^ 2> (fx) for different schemes, A^jg = 300 MeV and common normalization 



4 1/2) (m c ) = l. 



Figure [12| : B^ 2 ^ (/x) for different schemes, A^jg = 300 MeV and common normalization 



Bi 1/2) (m c ) = l. 

,(1/2) 



Figure [T3|: Bf' J (fi) for different schemes, A^ = 300 MeV and common normalization 



Bi 1/2 \m c ) = l. 

,(1/2 



Figure [14] : B\ (fi) for different schemes, A^ = 300 MeV and common normalization 



Bi 1/2 \m c ) = 1. 



Figure [T5| : 2 (/i) as a function of -Bg 1//2 ^(m c ) in the HV scheme for Aj^g = 300 MeV. 

88 



Figure [T6| : The ranges of ImA t as a function of m t for parameter ranges I (left) and II 
(right) given in appendix C. 



Figure 17| : The most important contributions to P* 1 / 2 ) and p( 3 / 2 ) as a function of /i for 
LO and the NDR scheme. m t , A^ and the various hadronic parameters Bi were 
taken at their central value given in appendix C. 



Figure [T8| : The ranges of e'/e in the NDR scheme as a function of A-^ for m t = 130 GeV 
and parameter ranges I (light grey) and II (dark grey) given in appendix C. The three 
pairs of e'/e plots correspond to hadronic parameters (a) B^ 2 \m c ) = B^ 2 \m c ) = 
1.5, (b) B ( i /2 \m c ) = 1.5, B<i /2 \m c ) = 1 and (c) ^ 1/2) (m c ) = Bf /2 \m c ) = 1, 
respectively. 



Figure 19: Same as fig. IS but for m t = 160 GeV. 



Figure |20| : Same as fig. |T| but for m t = 200 GeV. 
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